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1 Introduction 



1.1 Overview 

It is an exciting time now in high-energy physics. For many years, ever since the Standard Model 
was established in the late 1970's, the next logical question in the search for the basic laws of physics 
has been that of the mechanism by which the weak interaction gauge symmetry is spontaneously 
broken. This seemed at the time the one important gap that kept the Standard Model from being a 
complete theory of the strong, weak, and electromagnetic interactions |ll2l3j . Thirty years later, after 
many precision experiments at high-energy e^e~ and hadron colliders, this is still our situation. In 
the meantime, another important puzzle has been recognized, the fact that 80% of the mass in the 
universe is composed of 'dark matter', a particle species not included in the Standard Model. Both 
problems are likely to be solved by new fundamental interactions operating in the energy range of a 
few hundred GeV. Up to now, there is no evidence from particle physics for such new interactions. 
But, in the next few years, this situation should change dramatically. Beginning in 2008, the CERN 
Large Hadron CoUider (LHC) should give us access to physics at energies well above 1 TeV and thus 
should probe the energy region responsible for electroweak symmetry breaking. Over a longer term, 
we can look forward to precision experiments in e^e~ annihilation in this same energy region at the 
proposed International Linear Collider (ILC). 

Given this expectation, it is important for all students of elementary particle physics to form 
concrete ideas of what new phenomena we might find as we explore this new energy region. Of 
course, we have no way of knowing exactly what we will find there. But this makes it all the more 
important to study the alternative theories that have been put forward and to understand their 
problems and virtues. 

Many different models of new physics relevant to electroweak symmetry breaking are being 
discussed at this TASI school. Among these, supersymmetry has pride of place. Supersymmetry (or 
SUSY) provides an explicit realization of all of the aspects of new physics expected in the hundred 
GeV energy region. Because SUSY requires only weak interactions to build a realistic theory, it 
is possible in a model with SUSY to carry out explicit calculations and find the answers that the 
model gives to all relevant phenomenological questions. 

In these lectures, I will give an introduction to supersymmetry as a context for building models 
of new physics associated with electroweak symmetry breaking. Here is an outline of the material: 
In Section 2, I will develop appropriate notation and then construct supersymmetric Lagrangians 
for scalar, spinor, and vector fields. In Section 3, I will define the canonical phenomenological 
model of supersymmetry, the Minimal Supersymmetric Standard Model (MSSM). I will discuss the 
quantum numbers of new particles in the MSSM and the connection of the MSSM to the idea of 
grand unification. 

The remaining sections of these lectures will map out the phenomenology of the new particles 
and interactions expected in models of supersymmetry. I caution you that I will draw only those 
parts of the map that cover the simplest and most well-studied class of models. Supersymmetry 
has an enormous parameter space which contains many different scenarios for particle physics, more 
than I have room to cover here. I will at least try to indicate the possible branches in the path and 
give references that will help you follow some of the alternative routes. 

With this restriction, the remaining sections will proceed as follows: In Section 4, I will compute 
the mass spectrum of the MSSM from its parameters. I will also discuss the parameters of the 
MSSM that characterize supersymmetry breaking. In Section 5, I will describe how the MSSM 
parameters will be measured at the LHC and the ILC. Finally, Section 6 will discuss the answers 
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that supersymmetry gives to the two major questions posed at the beginning of this discussion, the 
origin of electroweak symmetry breaking, and the origin of cosmic dark matter. 

Although I hope that these lectures will be useful to students in studying supersymmetry, there 
are many other excellent treatments of the subject available. A highly recommended introduction 
to SUSY is the 'Supersymmetry Primer' by Steve Martin [B]. An excellent presentation of the 
formalism of supersymmetry is given in the texbook of Wess and Bagger [7]. Supersymmetry has 
been reviewed at previous TASI schools by Bagger [8], Lykken [9], and Kane [10], among others. Very 
recently, three textbooks of phenomenological supersymmetry have appeared, by Drees, Godbole, 
and Roy [TT], Binetruy \TZ, and Baer and Tata piJ:. A fourth textbook, by Dreiner, Haber, and 
Martin [Tl], is expected soon. 

It would be wonderful if all of these articles and books used the same conventions, but that is too 
much to expect. In these lectures, I will use my own, somewhat ideosyncratic conventions. These 
are explained in Section 2.1. Where possible, within the philosophy of that section, I have chosen 
conventions that agree with those of Martin's primer [B]. 



1.2 Motivation and Structure of Supersymmetry 

If we propose supersymmetry as a model of electroweak symmetry breaking, we might begin by 
asking: What is the problem of electroweak symmetry breaking, and what are the alternatives for 
solving it? 

Electroweak symmetry is spontaneously broken in the minimal form of the Standard Model, which 
I will refer to as the MSM. However, the explanation that the MSM gives for this phenomenon is 
not satisfactory. The sole source of symmetry breaking is a single elementary Higgs boson field. All 
mass of quarks, leptons, and gauge bosons arise from the couplings of those particles to the Higgs 
field. 

To generate symmetry breaking, we postulate a potential for the Higgs field 

V ^ ^^'\^\' + x\v\\ (1) 

shown in Fig. [l] The assumption that < is the complete explanation for electroweak symmetry 
breaking in the MSM. Since /z is a renormalizable coupling of this theory, the value of /i cannot be 
computed from first principles, and even its sign cannot be predicted. 

In fact, this explanation has an even worse feature. The parameter /i^ receives large additive 
radiative corrections from loop diagrams. For example, the two diagrams shown in Fig. [2] are 
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Figure 2: Two Standard Model diagrams that give divergent corrections to the Higgs mass parameter fi^. 



ultraviolet divergent. Supplying a momentum cutoff A, the two diagrams contribute 

If we view the MSM as an effective theory, A should be taken to be the largest momentum scale at 
which this theory is still valid. The presence of large additive corrections implies that the criterion 
/i^ < is not a simple condition on the underlying parameters of the effective theory. The radiative 
corrections can easily change the sign of fj?. Further, if we insist that the MSM has a large range 
of validity, the corrections become much larger than the desired result. To obtain the Higgs field 
vacuum expectation value required for the weak interactions, should be about 100 GeV. If we 
insist at the same time that the MSM is valid up to the Planck scale, A ~ 10^^ GeV, the formula 
([2]) requires a cancellation between the bare value of /i and the radiative corrections in the first 
36 decimal places. This problem has its own name, the 'gauge hierarchy problem'. But, to my 
mind, the absence of a logical explantion for electroweak symmetry breaking in the MSM is already 
problem enough. 

How could we solve this problem? There are two different strategies. One is to look for new 
strong-couplings dynamics at an energy scale of 1 TeV or below. Then the Higgs field could be 
composite and its potential could be the result, for example, of pair condensation of fermion con- 
stituents. Higgs actually proposed that his field was a phenomenological description of a fermion 
pair condensation mechanism similar to that in superconductivity Sometime later, Susskind ^ 
and Weinberg [3 proposed an explicit model of electroweak symmetry breaking by new strong in- 
teractions, called 'technicolor'. 

Today, this approach is disfavored. Technicolor typically leads to fiavor-changing neutral cur- 
rents at an observable level, and also typically confiicts with the accurate agreement of precision 
electroweak theory with experiment. Specific models do evade these difficulties, but they are highly 
constrained [S]. 

The alternative is to postulate that the electroweak symmetry is broken by a weakly-coupled 
Higgs field, but that this field is part of a model in which the Higgs potential is computable. In 
particular, the Higgs mass term /i^|<y3p should be generated by well-defined physics within the model. 
A prerequisite for this is that the term not receive additive radiative corrections. This requires 
that, at high energy, the appearance of a nonzero /i^ in the Lagrangian should be forbidden by a 
symmetry of the theory. 

There are three ways to arrange a symmetry that forbids the term /i^j^jp. We can postulate a 
symmetry that shifts f 

Sip — ev . (3) 

We can postulate a symmetry that connects ip to a gauge field, whose mass can then forbidden by 
gauge symmetry 

Sifi^e-A . (4) 
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We can postulate a symmetry that connects to a fermion field, whose mass can then be forbidden 
by a chiral symmetry. 

5lp — € ■ -ij} . (5) 

The options ([S]) and Q lead, respectively, to 'little Higgs' models |15I16I17] and to models with 
extra space dimensions ! 18119] , The third option leads to supersymmetry. This is the route we will 
now follow. 

The symmetry ([s]) looks quite innocent, but it is not. In quantum theory, a symmetry that links 
a boson with a fermion is generated by a conserved charge Qa that carries spin-1/2 

[Qc.,'A^^c, [Qa,i/] = 0. (6) 

Such a Qa implies the existence of a conserved 4- vector charge R„i defined by 

{Q„,Q^}-27™^i?„ (7) 

(It may not be obvious to you that there is no Lorentz scalar component in this anticommutator, 
but I will show this in Section 2.1.) The charge Rm is conserved, because both Q and commute 
with H . It is nonzero, as we can see by taking the expectation value of ([t]) in any state and setting 

a = l3 

{A\ {g., Ql} \A) = {A\ Q^Ql \A) + {A\ Q^qI \A) 

^\\Q^\A)\\' + \\Ql\A)\\' . (8) 

This expression is non-negative; it can be zero only if Qa and Qj^, annihilate every state in the 
theory. 

However, in a relativistic quantum field theory, we do not have the freedom to introduce arbitrary 
charges that have nontrivial Lorentz transformation properties. Conservation of energy-momentum 
and angular momentum are already very constraining. For example, in two-body scattering, the 
scattering amplitude for fixed center of mass energy can only be a function of one variable, the 
center of mass scattering angle 9. If one adds a second conserved 4- vector charge, almost all values 
of 9 will also be forbidden. Coleman and Mandula proved a stronger version of this statement: In 
a theory with an addtional conserved 4-vector charge, there can be no scattering at all, and so the 
theory is trivial [20 . 

If we would like to have ([s]) as an exact symmetry, then, the only possibility is to set Rm = Pm- 
That is, the square of the fermionic charge Qa must be the total energy-momentum of everything. 
We started out trying to build a theory in which the fermionic charge acted only on the Higgs field. 
But now, it seems, the fermionic charge must act on every field in the theory. Everything — quarks, 
leptons, gauge bosons, even gravitons — must have partners under the symmetry generated by Qa- 
Qa is fermionic and carries spin ^. Then every particle in the theory must have a partner with the 
opposite statistics and spin differing by ^ unit. 

The idea that the transformation (|5| leads to a profound generalization of space-time symmetry 
was discovered independently several times in the early 1970's |22l23j . The 1974 paper by Wess and 
Zumino |24j which gave simple linear realizations of this algebra on multiplets of fields launched the 
detailed exploration of this symmetry and its application to particle physics. 

The pursuit of ^ then necessarily leads us to introduce a very large number of new particles. 
This seems quite daunting. It might be a reason to choose one of the other paths, except that these 
also lead to new physics models of similarly high complexity. I encourage you to carry on with this 
line of analysis a bit longer. It will lead to a beautiful structure with many interesting implications 
for the theory of Nature. 
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2 Formalism of Supersymmetry 



2.1 Fermions in 4 Dimensions 

To work out the full consequences of ([s]), we will need to write this equation more precisely. 
To do this, we need to set up a formalism that describes relativistic fermions in four dimensions 
in the most general way. There is no general agreement on the best conventions to use, but every 
discussion of supersymmetry leans heavily on the particular choices made. I will give my choice of 
conventions in this section. 

There are two basic spin-| representations of the Lorentz group. Each is two-dimensional. The 
transformation laws are those of left- and right-handed Weyl (2-component) fermions, 

V'L->(l-ia-a/2-/3-a/2)^/.L 

V-fl^ (l-za-a/2 + /3-(7/2)^fl, , (9) 

where a is an infinitesimal rotation angle and f3 is an infinitesimal boost. The four-component 
spinor built from these ingredients, \l/ = (i/'Lj V'-r)) is a Dirac fermion. 



Define the matrix 



c=-^a^=r . (10) 



This useful matrix satisfies = —1, = —c. The combination 

1plLC'ip2L = -ea/3V'lLaV'2L/3 (H) 

is the basic Lorentz invariant product of spinors. Many treatments of supersymmetry, for example, 
that in Wess and Bagger's book [7], represent c implicitly by raising and lowering of spinor indices. 
I will stick to this more prosaic approach. 

Using the identity ac = —c{a)'^, it is easy to show that the quantity {—ctpD transforms like tpf^. 
So if we wish, we can replace every ipn by a V'l and write all fermions in the theory as left-handed 
Weyl fermions. With this notation, for example, we would call and fermions and and 
antifermions. This convention does not respect parity, but parity is not a symmetry of the Standard 
Model. The convention of representing all fermions in terms of left-handed Weyl fermions turns out 
to be very useful for not only for supersymmetry but also for other theories of physics beyond the 
Standard Model. 

Applying this convention, a Dirac fermion takes the form 



7'" -(in ) (13) 



Write the Dirac matrices in terms of 2 x 2 matrices as 


^ a"' 

with 

a™ = (1, ct)" ct™ = (1, -ct)™ cct™ = (ct")'^c (14) 
Then the Dirac Lagrangian can be rewritten in the form 

= ipli^iW ■ di/jiL + iplLia ■ dip2L 

-{rmljJi^cip2L - m*^pl^cip2L) ■ (15) 
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For the bilincars in the last hne, we can use fermion anticommutation and the antisymmetry of c 
to show 

V'ilCi/'2L = +'4)li^ci^iL ■ (16) 

and, similarly, 

{^^L0i^2L)^ = i^lLi-cWlL = -^1lCV'2L • (17) 

The mass term looks odd, because it is fermion number violating. However, the definition of fermion 
number is that given in the previous paragraph. The fields ifiiL and ■02L annihilate, respectively, e£ 
and ej. So this mass term generates the conversion of to e^, which is precisely what we would 
expect a mass term to do. 

If we write all fermions as left-handed Weyl fermions, the possibilities for fermion actions are 
highly restricted. The most general Lorentz-invariant free field Lagrangian takes the form 

L = ■i/'fc«CT • dipk - ^{mjki/jjctpk - fUj^ipyipl) . (18) 

where j, k index the fermion fields. Here and in the rest of these lectures, I drop the subscript L. 
The matrix rrijk is a complex symmetric matrix. For a Dirac fermion, 



as we have seen in (15). This matrix respects the charge 

Q^i = , QV2 = -^2 , (20) 

which is equivalent to the original Dirac fermion number. A Majorana fermion is described in the 
same formalism by the mass matrix 

rrijk = mSjk ■ (21) 

The most general fermion mass is a mixture of Dirac and Majorana terms. We will meet such 
fermion masses in our study of supersymmetry. These more general mass matrices also occur in 
other new physics models and in models of the masses of neutrinos. 

The SUSY charges are four-dimensional fermions. The minimum set of SUSY charges thus 
includes one Weyl fermion Qa and its Hermitian conjugate Q^^. We can now analyze the anticom- 
mutator {QcQ]^}- Since the indices belong to different Lorentz representations, this object does 
not contain a scalar. The indices transform as do the spinor indices of ct™, and so we can rewrite 
^ with i?™ = P" as 

{Q„,g^} = 2<^P,„ . (22) 

It is possible to construct quantum field theories with larger supersymmetry algebras. These 
must include (22 1, and so the general form is pi] 

{Ql,Qy}^2aZ'pPm6'' , (23) 
for i, j = 1 . . . N . This relation can be supplemented by a nontrivial anticommutator 

{Q]^,Q^^}^2e^p&^ (24) 

where the central charge Q'^ is antisymmetric in [ij]. Theories with > 4 necessarily contain 
particles of spin greater than 1. Yang- Mills theory with N — 4 supersymmetry is an especially 
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beautiful model with exact scale invariance and many other attractive formal properties [25j . In 
these lectures, however, I will restrict myself to the minimal case of = 1 supersymmetry. 

I will discuss supersymmetry transformations using the operation on fields 

S^^=[ecQ + Q^c^*,<^>] . (25) 

Note that the operator contains pairs of anticommuting objects and so obeys commutation rather 
than anticommutation relations. The operator P„j acts on fields as the generator of translations, 
Pjn = idjn- Using this, we can rewrite ( |22| ) as 

[<5^,<5,] = 2z(cV'"77-77V'"^) dm (26) 

I will take this equation as the basic (anti)-commutation relation of supersymmetry. In the next two 
sections, I will construct some representations of this commutation relation on multiplets of fields. 



2.2 Supersymmetric Lagrangians with Scalars and Fermions 



The simplest representation of the supersymmetry algebra ( 26 1 directly generalizes the transfor- 
mation ([5]) from which we derived the idea of supersymmetry. The full set of fields required includes 
a complex- valued boson field (p and a Weyl fermion field ip. These fields create and destroy a scalar 
particle and its antiparticle, a left-handed massless fermion, and its right-handed antiparticle. Note 
that the particle content has an equal number of fermions and bosons. This particle content is called 
a chiral supermultiplet. 

I will now write out the transformation laws for the fields corresponding to a chiral supermultiplet. 
It is convenient to add a second complex- valued boson field F that will have no associated particles. 
Such a field is called an auxiliary field. We can then write the transformations that generalize (|5| as 

S^(j) = %/2^'^ctA 

5^i; = V2ic7''cCdn(t) + V2F^ 

S^F = -V2i£}a'^dmi^ . (27) 
The conjugates of these transformations are 

S^cf)* = -\/2V'^cr 

S^F* = V2id„,tp^a"'^ . (28) 
These latter transformations define the antichiral supermultiplet. I claim that the transformations 



(27 1 and (28 1, first, satisfy the fundamental commutation relation (26) and, second, leave a suitable 
Lagrangian invariant. Both properties are necessary, and both must be checked, in order for a set 
of transformations to generate a symmetry group of a field theory. 



The transformation laws (|27| seem complicated. You might wonder if there is a formalism that 
generates these relations automatically and manipulates them more easily than working with the 
three distinct component fields {(f>,4^,F). In the next section, I will introduce a formalism called 
superspace that makes it almost automatic to work with the chiral supermultiplet. However, the 
superspace description of the multiplet containing gauge fields is more complicated, and the difficulty 
of working with superspace becomes exponentially greater in theories that include gravity, higher 
dimensions, or iV > 1 supersymmetry. At some stage, one must go back to components. I strongly 
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recommend that you gain experience by working through the component field calculations described 
in these notes in full detail, however many large pieces of paper that might require. 



To verify each of the two claims I have made for (27 1 requires a little calculation. Here is the 
check of the commutation relation applied to the field 0: 

= 2z[^V>-77t^"^]9„^ 

(29) 

The check of the commutation relation applied to F is equally straightforward. The check on ?/; is a 
bit lengthier. It requires a Fierz identity, that is, a spinor index rearrangment identity. Specifically, 
we need 



(30) 



which you can derive by writing out the four components explicitly. After some algebra that involves 
the use of this identity, you can see that the SUSY commutation relation applied to ip also takes the 
correct form. 



Next, I claim that the Lagrangian 



F*F 



(31) 



is invariant to the transformation (27 1. I will assume that the Lagrangian (31 1 is integrated J d^x 
and use integration by parts freely. Then 

JjL = 9"V*a™(%/2C^cV') + (-V2a"V^cr)50 
+'4)ha ■ d[V2icr''cCd^(l> + V2^F] 
+ [y/2idn4>*(^ca'' + y/2S} F*]i<j ■ dip 
+F*[~V2i£}a"'d,nip] + [V2id^iP^a"'i]F 
= -0*v^e^c9V + y25„(^*e^ca"(f™a„a„^ 
+^f2ip'^cCd^<t) - \/2V'V'"fT"cr9™9„</) 
+V2ii)^a"'dmFi + V2idm'>p'^a"'£,F 
-V2iS}F*a"'draip + V2iF*^^a"'dmip 



= 



(32) 



In the final expression, the four lines cancel line by line. In the first two lines, the cancellation is 
made by using the identity {a ■ d){a ■ d) = d^. 

So far, our supersymmetry Lagrangian is just a massless free field theory. However, it is possible 
to add rather general interactions that respect the symmetry. Let W{4') be an analytic function of 
(j), that is, a function that depends on (f> but not on (f>* . Let 



F 



dW 



1 ,7, B^W 
2^ '"^^ 



(33) 



I claim that L^i/ is invariant to (27 1. Then we can add {hw + ^^^^ field Lagrangian to 

introduce interactions into the theory. The function W is called the superpotential. 
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We can readily check that hw is indeed invariant: 



-Vctp—j^V^fc^ . (34) 



The second hne rearranges to 



(35) 



which is a total derivative. The third line is proportional to ipa'^'fj'ip-y, which vanishes by fermion 
antisymmetry since the spinor indices take only two values. Thus it is true that 

6^Lw = . (36) 



The proofs of invariance that I have just given generalize straightforwardly to systems of several 
chiral supermultiplets. The requirement on the superpotential is that it should be an analytic 
function of the complex scalar fields (j)k- Then the following Lagrangian is supersymmetric: 



d'^^ldrnA + ^lio ■ d^k + F*Fk + Lw 



t 



where 



JW 



dW 



1 y 



' d4>k ^ " ocpjocpk 

In this Lagrangian, the fields are Lagrange multipliers. They obey the constraint equations 

dW 



Ft 



(37) 
(38) 

(39) 



Using these equations to eliminate the Fk, we find an interacting theory with the fields (pk and V^fc, a 
Yukawa coupling term proportional to the second derivative of W , as given in (38 1, and the potential 
energy 



dW 



(40) 



I will refer to Vp as the F-term potential. Later we will meet a second contribution Vd, the D-term 
potential. These two terms, both obtained by integrating out auxiliary fields, make up the classical 
potential energy of a general supersymmetric field theory of scalar, fermion, and gauge fields. 

The simplest example of the F-term potential appears in the theory with one chiral supermultiplet 
and the superpotential W — ^m<j)'^. The constraint equation for F is [27] 



F* — ~m 

After eliminating F, we find the Lagrangian 
L 



(41) 



(42) 
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This is a theory of two free scalar bosons of mass |m| and a free Majorana fermion with the same 
mass \m\. The Majorana fermion has two spin states, so the number of boson and fermion physical 
states is equal, as required. 



The form of the expression (40 1 implies that Vp > 0, and that Vp = only if all Fk = 0. This 



constraint on the potential energy follows from a deeper consideration about supersymmetry. Go 



back to the anticommutation relation (22), evaluate it for a = f3, and take the vacuum expectation 
value. This gives 

(0| {Qa, |0) = (0| {H - P3) |o) = (0| H |0) , (43) 

since the vacuum expectation value of vanishes by rotational invariance. Below ([7|, I argued 
that the left-hand side of this equation is greater than or equal to zero. It is equal to zero if and 
only if 

Q„ |0) = Qi |0) - (44) 

The formulae (44 1 give the criterion than the vacuum is invariant under supersymmetry. If this 
relation is not obeyed, supersymmetry is spontaneously broken. Taking the vacuum expectation 
value of the transformation law for the chiral representation, we find 

(0| [^cQ + Q^cf , Vfe] |0) = (0| V2za"ra„0fc + ^ |0) 

= ^(0|F,.|0) . (45) 

In the last line I have used the fact that the vacuum expectation value of (/>(x) is translation invariant. 



so its derivative vanishes. The left-hand side of (45 1 vanishes if the vacuum state is invariant under 
supersymmetry. 

The results of the previous paragraph can be summarized in the following way: If supersymmetry 
is a manifest symmetry of a quantum field theory, 

(0| H\0) = , and (0| Fk |0) = (46) 

for every F field of a chiral multiplet. In complete generality, 

(0|ff|0)>0. (47) 

The case where (H) is positive and nonzero corresponds to spontaneously broken supersymmetry. If 



the theory has a state satisfying (44 1, this is necesssarily the state in the theory with lowest energy. 
Thus, supersymmetry can be spontaneously broken only if a supersymmetric vacuum state does not 
exisQ 

For the moment, we will work with theories that preserve supersymmetry. I will give examples 
of theories with spontaneous supersymmetry breaking in Section 3.5. 

The results we have just derived are exact consequences of the commutation relations of super- 
symmetry. It must then be true that the vacuum energy of a supersymmetric theory must vanish 



in perturbation theory. This is already nontrivial for the free theory (42 1. But it is correct. The 
positive zero point energy of the boson field exactly cancels the negative zero point energy of the 
fermion field. With some effort, one can show the cancellation also for the leading-order diagrams 
in an interacting theory. Zumino proved that this cancellation is completely general [29j. 

I would like to show you another type of cancellation that is also seen in perturbation theory in 
models with chiral fields. Consider the model with one chiral field and superpotential 

W^^ct)^ . (48) 



*It is possible that a supersymmetric vacuum state might exist but that a higher-energy vacuum state might be 
metastable. A model built on this metastable state would show spontaneous breaking of supersymmetry |26| . 
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b.) 




Figure 3: Perturbation theory for the supersymmetric model ( |49[ ): (a) vertices of the model; (b) corrections 
to the fermion and scalar masses. 



\ 
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F + * 




(|) + X)) 





= 



Figure 4: Scheme of cancellations of one-loop corrections to the F-term potential. 



After eliminating F, the Lagrangian becomes 

L = a0*<9,„0 + V^iCT-9V-A(0V'^c7/'-</)*V^cV'*)-A2|0|4 . (49) 

The vertices of this theory are shown in Fig.|3|a). 

From our experience in (|2|, we might expect to find an addditive radiative correction to the 
scalar mass. The corrections to the fermion and scalar mass terms are given by the diagrams in 
Fig.|3|b). Actually, there are no diagrams that correct the fermion mass; you can check that there it 
is not possible to match the arrows appropriately. For the scalar mass correction, the two diagrams 
shown contribute 



la ■ p la 



pi 



(50) 



Using a -pa-p = p^ in the second term and then taking the trace, we see that these two contributions 
cancel precisely. In this way, supersymmetry really does control radiative corrections to the Higgs 
mass, following the logic that we presented in Section 1.2. 

In fact, it can be shown quite generally that not only the mass term but the whole superpotential 
W receives no additive radiative corrections in any order of perturbation theory 30J. For example, 
the one-loop corrections to quartic terms in the Lagrangian cancel in a simple way that is indicated 
in Fig. [4] The field strength renormalization of chiral fields can be nonzero, so the form of W can be 
changed by radiative corrections by the rescaling of fields. Examples are known in which W receives 
additive radiative corrections from nonperturbative effects |31j . 
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2.3 Superspace 



Because the commutation relations of supersymmetry include the generators of translations, 
supersymmetry is a space-time symmetry. It is an attractive idea that supersymmetry is the natural 
set of translations on a generalized space-time with commuting and anticommuting coordinates. In 
this section, I will introduce the appropriate generalization of space-time and use it to re-derive some 
of the results of Section 2.2. 

Consider, then, a space with four ordinary space-time coordinates x'' and four anticommuting 
coordinates {Oa,Oa)- I will take the coordinates 9a to transform as 2-component Weyl spinors; the 
9a are the complex conjugates of the 9^- This is superspace. A superfield is a function of these 
superspace coordinates: ^{x,9,9). 

It is tempting to define supersymmetry transformations as translations 9^9 + ^. However, this 
does not work. These transformations commute, [<5{,<577] = 0, and we have seen in Section 1.2 that 
this implies that the S-matrix of the resulting field theory must be trivial. To construct a set of 
transformations with the correct commutation relations, we must write 

S^<^ = , (51) 

where 

Q« = - i9a"'d„?j ^ + ia^9d„?j . (52) 

This is a translation of the fermionic coordinates {9, 9) plus a translation of the ordinary space-time 
coordinates proportional to 6*, 9. It is straightforward to show that these operators satisfy 

[Q^, Q„] = -2^ (CV'"77 - j/t^™^) . (53) 
Despite the fact that this equation has an extra minus sign on the right-hand side with respect to 



(26 1, it is the relation that we want. (The difference is similar to that between active and passive 



transformations.) Combined with the decomposition of the superfield that I will introduce below. 



this relation will allow us to derive the chiral supermultiplet transformation laws ( 27 ) 



Toward this goal, we need one more ingredient. Define the superspace derivatives 

D^=^- i{9o"')o.dm Da = -^ + i{a"'9)ad^ , (54) 

Cfa 09a 

such that {Da^y = Da^^- These operators commute with Q^: 

[Da,Qi] = [Da,Qi] = 0. (55) 

Thus, we can constrain <i> by the equation 

= or Da<^ = , (56) 

and these constraints are consistent with supersymmetry. What we have just shown is that the 
general superfield ^{x,9,9) is a reducible representation of supersymmetry. It can be decomposed 



into a direct sum of three smaller representations, one constrained by the first of the relations (56 1, 
one constrained by the second of these relations, and the third containing whatever is left over in $ 
when these pieces are removed. 

Let's begin with the constraint = 0. The solution of this equation can be written 

$(x, 9, 9) = <^{x + i9a"'9, 9) , (57) 
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that is, this solution is parametrized by a general function of x and 6. Since is a two-component 
anticommuting object, this general function of x and 6 can be represented as 



0) = (l){x) + V2e^ci^{x) + 9^c9F{x) . 



(58) 



The field content of this expression is exactly that of the chiral supermultiplet. The supersymmetry 
transformation of this field should be 



= Q^<^{x + iOa''- 



(59) 



It is straightforward to compute the right-hand side of (59 1 in terms of 0, 0, and the component 



fields of ( 58 1 . The coefficients of powers of are precisely the supersymmetry variations given in 



( 27 1 . Thus a superfield satisfying 



(60) 



is equivalent to a chiral supermultiplet, and the transformation (59 1 gives the supersymmetry trans- 



formation of this multiplet. A superfield satisfying (60 1 is called a chiral superfield. Similarly, a 
superfield satisfying 

= (61) 

is called an antichiral superfield. This superfield has a component field decomposition ((jj* ^tp* , F*), 
on which induces the transformation ( 28 1 . I will describe the remaining content of the general 
superfield (f> in Section 2.5. 

A Lagrangian on Minkowski space is integrated over d'^x. A superspace Lagrangian should be 
also be integrated over the 6 coordinates. Integration over fermionic coordinates is defined to be 
proportional to the coefficient of the highest power of 9. I will define integration over superspace 
coordinates by the formulae 



d-'ei 







j d^e{9'^c6 



(62) 



and their conjugates. To use these formulae, expand the superfields in powers of 6 and pick out the 
terms proportional to {9^ 
analytic function of <i>. 



Then, if $ is a chiral superfield constrained by (60) and T4^(<i>) is an 



d'^9 Wi<i>) = F{x) 



dW 



1 ,r ,d'^W 
2^ 



(63) 



where, in the second line, W on the right-hand side is evaluated with <t> = 0(a;). With somewhat 
more effort, one can show 



d^e / d^9 = + i)\u ■ d^p + f*f 



(64) 



These formulae produce the invariant Lagrangians of chiral supermultiplets from a superspace 
point of view. The most general Lagrangian of chiral superfields takes the form 



d'^9K{<i>,^^) 



d^9W{^) + / d^9{W{^)) 



(65) 



where VK(<i>) is an analytic function of complex superfields and K{^,^^) is a general real-valued 
function of the superfields. The Lagrangian (37 1 is generated from this expression by taking 
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ii:($,$t) ^ ^l^k- The most general renormalizahle Lagrangian of chiral supermultiplets is ob- 
tained by taking K to be of this simple form and taking to be a polynomial of degree at most 3. 



Because the integral (PO exposes the Lagrange multiplier in (58 1, I will refer to a term with 
this superspace integral as an F-term. For similar reasons that will become concrete in the next 
section, I will call a term with a d'^O integral a D-term. 

In the remainder of these lectures, I will restrict myself to discussing renormalizable supersym- 
metric theories. But, still, it is interesting to ask what theories we obtain when we take more general 
forms for K. The Lagrangian for (j) turns out to be a nonlinear sigma model for which the target 
space is a complex manifold with the metric [32] 

A complex manifold whose metric is derived from a potential in this way is called a Kdhler manifold. 
The function K is the Kdhler potential. It is remarkable that, wherever in ordinary quantum field 
theory we find a general structure from real analysis, the supersymmetric version of the theory has 
a corresponding complex analytic structure. 

Now that we have a Lagrangian in superspace, it is possible to derive Feynman rules and compute 
Feynman diagrams in superspace. I do not have space here to discuss this formalism; it is discussed, 
for example, in and [30 . I would like to state one important consequence of this formalism. It 
turns out that, barring some special circumstances related to perturbation theory anomalies, these 
Feynman diagrams always generate corrections to the effective Lagrangian that are D-terms, 



d''9X{^,<^^) . (67) 

The perturbation theory does not produce terms that are integrals J d^d. This leads to an elegant 
proof of the result cited at the end of the previous section that the superpotential is not renormalized 
at any order in perturbation theory j30j. 



2.4 Supersymmetric Lagrangians with Vector Fields 

To construct a supersymmetric model that can include the Standard Model, we need to be able 
to write supersymmetric Lagrangians that include Yang-Mills vector fields. In this section, I will 
discuss how to do that. 

To prepare for this discussion, let me present my notation for gauge fields in a general quantum 
field theory. The couplings of gauge bosons to matter are based on the covariant derivative, which 
I will write as 

V„A = (9„, - igAlt-^)(b (68) 

for a field </> that belongs to the representation R of the gauge group G. In this formula, tfj are the 
representation matrices of the generators of G in the representation R. These obey 

[fkAJ-^r'^tR (69) 

The coefficients f""^^ are the structure constants of G. They are independent of R; essentially, their 
values define the multiplication laws of G. They can be taken to be totally antisymmetric. 

The generators of G transform under G according to a representation called the adjoint repre- 
sentation. I will denote this representation by i? = G. Its representation matrices are 

= z/"'^^ (70) 
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These matrices satisfy ( 69 ) by virtue of the Jacobi identity. The covariant derivative on a field in 



the adjoint representation takes the form 

I?,„<i>" = a™*" + gr'"'Al<^>^ (71) 
The field strengths are defined from the covariant derivative (in any representation) by 

[V„M^~^9KJr- (72) 

This gives the familiar expression 

F^^ = ~ d^Al^ + gr'^AlAl . (73) 

Now we would like to construct a supersymmetry multiplet that contains the gauge field A"^. 
The fermion in the multiplet should differ in spin by ^ unit. To write a renormalizable theory, we 
must take this to be a spin-^ Weyl fermion. I will then define the vector supermultiplet 

{Al^Al.D^) (74) 

including the gauge field, a Weyl fermion in the adjoint representation of the gauge group, and an 
auxililary real scalar field, also in the adjoint representation, that will have no independent particle 
content. The particle content of this multiplet is one massless vector boson, with two transverse 
polarization states, and one massless fermion and antifermion, for each generator of the gauge group. 
The fermion is often called a gaugino. The number of physical states is again equal between bosons 
and fermions. 

The supersymmetry transformations for this multiplet are 

S^D'^ = -zl^V^P^A"^ - P^AtV™?] (75) 

where 

ct"" = |(cr™CT" - cr"CT™) . (76) 

I encourage you to verify that these tranformations obey the algebra 

[5^ , <S„] = 2i (e V"77 - 7? dm + <5a , (77) 
where 5a is a gauge tranformation with the gauge parameter 

a = -2z(eV",y-,7V"0^:^ • (78) 



Acting on A°, the extra term 5a in (|77| can be combined with the translation to produce the 
commutation relation 

[5^,5,,] - 2i (^t^™^ „ ^t^™^) (p^^xr . (79) 

This rearrangement applies also for the auxiliary field and for any matter field that tranforms 
linearly under G. The gauge field A°'"^ does not satisfy this last criterion; instead, we find 

[5^,5.,]Al, = 2z(^t^n^ _ r;ta"0(9nA"„ - 

^2i{£}a-r^-r,^a-£,)Kra (80) 
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The proof that ( 75 1 satisfies the supersymmetry algebra is more tedious than for ( 29 1 , but it is not 

ition 

(81) 



actually difficult. For the transformation of A" we need both the Fierz identity (30 1 and the relation 
The matrices cr^'c and ca^'^ are symmetric in their spinor indices. 



Again, the transformation laws leave a simple Lagrangian invariant. For the vector supermulti- 
plet, this Lagrangian is that of the renormalizable Yang-Mills theory including the gaugino: 



The kinetic term for Z?" contains no derivatives, so this field will be a Lagrange multiplier. 



(82) 



The vector supermultiplet can be coupled to matter particles in chiral supermultiplets. To do this, 
we must first modify the transformation laws of the chiral supermultiplet so that the commutators 



of supersymmetry transformations obey ( |77| or (79 1. The modified transformation laws are: 

S^ijj = V2ia"cCVn(t)+V2F^ 



(83) 



In this formula, the chiral fields 0, V'l F must belong to the same representation of G, with 1°" 
a representation matrix in that representation. From the transformation laws, we can construct 
the Lagrangian. Start from (31 1, replace the derivatives by covariant derivatives, add terms to the 
Lagrangian involving the A" to cancel the supersymmetry variation of these terms, and then add 
terms involving Z?" to cancel the remaining supersymmetry variation of the A° terms. The result is 



Ln =1?" 



(84) 



The proof that this Lagrangian is supersymmetric, 5jL = 0, is completely straightforward, but it 
requires a very large sheet of paper. 

The gauge invariance of the theory requires the superpotential Lagrangian Lvk to be invariant 
under G as a global symmetry. Under this condition, L^y, which contains no derivatives, is invariant 
under (83) without modification. The combination oihp, i^D, and Liy, with W a polynomial of 
degree at most 3, gives the most general renormalizable supersymmetric gauge theory. 

As we did with the F field of the chiral multiplet, it is interesting to eliminate the Lagrange 
multiplier D". For the Lagrangian which is the sum of (82) and (84), the equation of motion for D"" 
is 

D" = ~g(t)*t''(t) . (85) 

Eliminating gives a second potential energy term proportional to (D°)^. This is the D-term 
potential promised below (40). I will write the result for a theory with several chiral multiplets: 

\ 2 



(86) 



As with the F-term potential, Vd > and vanishes if and only if D° = 0. It can be shown by an 

(0|L»"|0)=0 (87) 



argument similar to ( 45 1 that 
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unless supersymmctry is spontaneously broken. 

It makes a nice illustration of this formalism to show how the Higgs mechanism works in super- 
symmetry. For definiteness, consider a supersymmetric gauge theory with the gauge group U{1). 

Introduce chiral supermultiplets and X, with charges +1, —1, and 0, respectively, and 

the superpotential 

W = A(0+(/)_ - v^)X . (88) 

The F — equations are 

F*x = (0+0- -v^) = Fl = ^±X = 0. (89) 
To solve these equations, set 

X = cb+ = v/y q^-^vy, (90) 

where y is a, complex- valued parameter. The D ~ equation is 

(t>l<f>+ ~ (f>U^ ^ . (91) 

This implies |?;| = 1. So y is a pure phase and can be removed by a t/(l) gauge transformation. 

Now look at the pieces of the Lagrangian that give mass to gauge bosons, fermions, and scalars. 
The gauge field receives mass from the Higgs mechanism. To compute the mass, we can look at the 
scalar kinetic terms 

0t(-2?')0+ + 0U-2?')0- - • • • + 4>Wa')^+ + 0-(.9'^')0- ■ (92) 
Putting in the vacuum expectation values 0+ = 4>- = "^c find 

= Ag^v^ (93) 
for the vector fields. The mode of the scalar field 

S(j)+ = r]/V2 5<j)^ = -77/\/2 , (94) 
with 77 real, receives a mass from the D-term potential energy 

^(40+-0t_^_)2 (95) 

Expanding to quadratic order in 77, we see that 77 also receives the mass m? = Ag^v^. The corre- 
sponding mode for 77 imaginary is the infinitesimal version of the phase rotation of y that we have 
already gauged away below (91 1. The mode of the fermion fields 

5ij+ = x/V2 <5^_ = -x/y2 (96) 

mixes with the gaugino through the term 

-\/2.g((?!)^A^C7/;+-0LA^cV'-) + /i.c. (97) 

Putting in the vacuum expectation values 0+ = 0- ^ v, we find a Dirac mass with the value 

TO = 2gv (98) 

In all, we find a massive vector boson, a massive real scalar, and a massive Dirac fermion, all with 
the mass to = 2gv. The system has four physical bosons and four physical fermions, all with the 
same mass, as supersymmetry requires. 
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2.5 The Vector Supermultiplet in Superspace 



The vector supermultiplet has a quite simple representation in superspace. This multiplet turns 
out to be the answer to the question that we posed in our discussion of superspace in the previous 
section: When the chiral and antichiral components of a general superfield are removed, what is left 
over? To analyze this issue, I will write a Lagrangian containing a local symmetry that allows us to 
gauge away the chiral and antichiral components of this superfield. Let V{x,9,6) be a real-valued 
superfield, acted on by a local gauge transformation in superspace 



(99) 



where A is a chiral superfield and A''' is its conjugate. Since A satisfies (60 1, its expansion in powers 
of 9 contains 

A(a;, 9, 9) = K{x + i9lf9, 9) = a{x) + ■■■ + i9a'''9d,na{x) + ■■■ (100) 
The general superfield V contains a terrrj^ 



V{x,9,9) 



+ 26l(f™6'A„(a;) + 



(101) 



So the superfield V contains a space-time vector field Am{x), and under (99 1, Am transforms as 



5 Am = -a,„(Rea) 
9 



(102) 



This is just what we would like for an Abelian gauge field. So we should accept (99 1 as the general- 
ization of the Abelian gauge transformation to superspace. 



The real- valued superfield transforming under (99) is called a vector superfield. To understand 



its structure, use the gauge transformation to remove all components with powers of 9 or 9 only. 
This choice is called Wess-Zumino gauge [33 . What remains after this gauge choice is 

V{x, 9, 9) = 29a"'9 Am{x) + 2f9^c\ - 29^9^ cX* + 9^fD . (103) 
This expression has exactly the field content of the Abelian vector supermultiplet {Am, A, -D). 

This gauge multiplet can be coupled to matter described by chiral superfields. For the moment, 
I will continue to discuss the Abelian gauge theory. For a chiral superfield <& with charge Q, the 
gauge transformation 

= iQA^ (104) 

contains a standard Abelian gauge transformation with gauge parameter Rea(x) and also preserves 
the chiral nature of <&. Then the superspace Lagrangian 



(105) 



is gauge-invariant. Using the representation (1031 and the rules (62 1, it is straightforward to carry 



out the integrals explicitly and show that ( 105 1 reduces to (84 1, with t°- — Q for this Abelian theory. 



We still need to construct the pure gauge part of the Lagrangian. To do this, first note that, 
because a quantity antisymmetrized on three Weyl fermion indices vanishes. 



D^D^X = 



(106) 



tThe factor 2 in this equation is convenient but disagrees with some standard treatments, e.g., [?]■ 
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2 

for any superfield X . Thus, acting with D makes any superfield a chiral superfield. The following 
is a chiral superfield that also has the property that its leading component is the gaugino field X(x): 

Wc. = - ^d\Dc)c.V . (107) 

Indeed, working this out in full detail, we find that Wa — Wa{x + iOaO, 6), with 

Wa{x, 0)^Xa + [(^a""F„,„ + D)9]^ + d'^ cO [a,„A*za"c]„ . (108) 

The chiral superfield Wa is the superspace analogue of the electromagnetic field strength. The 
Lagrangian 

d^e ]^W'^cW (109) 



reduces precisely to the Abelian version of ( 82 1 . It is odd that the kinetic term for gauge fields is an 
F-term rather than a D-term. It turns out that this term can be renormalized by loop corrections 
as a consequence of the trace anomaly [34J . However, the restricted form of the correction has 
implications, both some simple ones that I will discuss later in Section 4.3 and and more profound 
implications discussed, for example, in [35'36\ 

I will simply quote the generalizations of these results to the non-Abelian case. The gauge 
transformation of a chiral superfield in the representation R of the gauge group is 

$ ^ 6*^^°*°$ $t ^ $te-«At"f ^ (110) 

where A° is a chiral superfield in the adjoint representation of G and is is the representation of 
the generators of G in the representation R. The gauge transformation of the vector superfield is 

Then the Lagrangian 

' (f9(fe<i>''eS^'^^°^ (112) 



is locally gauge-invariant. Carrying out the integrals in the gauge (103) reduces this Lagrangian to 
([84]). 

The form of the field strength superfield is rather more complicated than in the Abelian case, 

WX = -5^S'e-s^°*°(i?c)„e9^°*° (113) 



In Wess-Zumino gauge, this formula does reduce to the non-Abelian version of (108), 

W^{x, 0) = AS -t- [(ia'»"F„"„, + D'^)e]^ + d^cd [I?„A-*a'"c]„ . (114) 
Then the Lagrangian 

d^dix\y/'^cW\ (115) 



reduces neatly to ( 82 ) . 

The most general renormalizable supersymmetric Lagrangian can be built out of these ingre- 
dients. We need to put together the Lagrangian (115), plus a term (112) for each matter chiral 
superfield, plus a superpotential Lagrangian to represent the scalar field potential energy. These 
formulae can be generalized to the case of a nonlinear sigma model on a Kahler manifold, with the 
gauge symmetry associated with an isometry of this target space. For the details, see [J- 
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2.6 R-Symmetry 



The structure of the general superspace action for a renormahzable theory of scalar and fermion 
fields suggests that this theory has a natural continuous symmetry. 

The superspace Lagrangian is 

L = J d^9tr[W'^cW]+ J d'*6'$^e»^*$ + J d^eW{<i>)+ J . (116) 

Consider first the case in which W{(j)) contains only dimensionless parameters and is therefore a 
cubic polynomial in the scalar fields. Then L is invariant under the U{1) symmetry 

<^k{x, 9) e-*2"/3$fc(x, e*"6') , Vix, 9, 9) ^{x, e'°'9, e"'"^) (117) 

or, in components, 

and the gauge fields are invariant. This transformation is called R- symmetry. Under R-symmetry, 
the charges of bosons and fermions differ by 1 unit, in such a way that that the gaugino and 
superpotential vertices have zero net charge. 



Since all left-handed fermions have the same charge under (1181, the R-symmetry will have an 
axial vector anomaly. It can be shown that the R-symmetry current (of dimension 3, spin 1) forms 
a supersymmetry multiplet together with the supersymmetry current (dimension |, spin |) and the 
energy-momentum tensor (dimension 4, spin 2) ^37j . All three currents have perturbation-theory 
anomalies; the anomaly of the energy-momentum tensor is the trace anomaly, associated with the 
breaking of scalar invariance by coupling constant renormalization. The R-current anomaly is thus 
connected to the running of coupling constants and gives a useful formal approach to study this 
effect in supersymmetric models. 



It is often possible to combine the transformation (117 1 with other apparent U{1) symmetries of 



the theory to define a non-anomalous U{1) R-symmetry. Under such a symmetry, we will have 

$fe(x, 9) e-^'^^-$fc(x, e^'e) , such that W{x, 9) ^ e^'^Wix, e"^e) . (119) 
Such symmetries also often arise in models in which the superpotential has dimensionful coefficients. 

In models with extended, > 1, supersymmetry, the R-symmetry group is also extended, to 
SU{2) for A = 2 and to SU{A) for A^ = 4 supersymmetry 



3 The Minimal Supersymmetric Standard Model 

3.1 Particle Content of the Model 

Now we have all of the ingredients to construct a supersymmetric generalization of the Standard 
Model. To begin, let us construct a version of the Standard Model with exact supersymmetry. To 
do this, we assign the vector fields in the Standard Model to vector supermultiplets and the matter 
fields of the Standard Model to chiral supermultiplets. 

The vector supermultiplets correspond to the generators of SU{3) x SU{2) x U{1). In these 
lectures, I will refer to the gauge bosons of these groups as AJ^, W^, and -B,„, respectively. I will 
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represent the Weyl fermion partners of these fields as g", w", fe. I will call these fields the gluino, 
wino, and bino, or, collectively, gauginos. In the later parts of these lectures, I will drop the tildes 
over the gaugino fields when they are not needed for clarity. 

I will assign the quarks and leptons to be fermions in chiral superfields. I will use the convention 
presented in Section 1.3 of considering left-handed Weyl fermions as the basic particles and right- 
handed Weyl fermions as their antiparticles. In the Standard Model, the left-handed fields in a 
fermion generation have the quantum numbers 

The field e is the left-handed positron; the fields u, d are the left-handed antiquarks. The right- 
handed Standard Model fermion fields are the conjugates of these fields. To make a generalization 



to supersymmetry, we will extend each of the fields in (1201 — for each of the three generations — to 
a chiral supermultiplet . I will use the symbols 

L ^ Q TL d (121) 

to represent both the supermultiplets and the scalar fields in these multiplets. Again, I will drop 
the tilde if it is unambiguous that I am referring to the scalar partner rather than the fermion. The 
scalar particles in these supermultiplets are called sleptons and squarks, collectively, sfermions. 

What about the Higgs field? The Higgs field of the Standard Model should be identified with 
a complex scalar component of a chiral supermultiplet. But it is ambiguous what the quantum 
numbers of this multiplet should be. In the Standard Model, the Higgs field is a color singlet with 
/ = |, but we can take the hypercharge of this field to be either Y — +^ or Y — depending on 
whether we take the positive hypercharge field or its conjugate to be primary. In a supersymmetric 
model, the choice matters. The superpotential is an analytic function of superfields, so it can only 
contain the field, not the conjugate. Then different Higgs couplings will be allowed depending on 
the choice that we make. 

The correct solution to this problem is to include both possibilities, That is, we include a Higgs 
supermultiplet with Y — + ^ and a second Higgs supermultiplet with ^ = — |- I will call the scalar 
components of these multiplets Hu and Hd, respectively: 



I will refer to the Weyl fermion components with these quantum numbers as h^- These fields or 
particles are called Higgsinos. 

I will argue below that it is necessary to include both Higgs fields in order to obtain all of 
the needed couplings in the superpotential. However, there is another argument. The axial vector 
anomaly of one U{1) and two SU{2) currents (Fig. [s]) must vanish to maintain the gauge invariance 
of the model. In the Standard Model, the anomaly cancels nontrivially between the quarks and 
the leptons. In the supersymmetric generalization of the Standard Model, each Higgsino makes a 
nonzero contribution to this anomaly. These contributions cancel if we include a pair of Higgsinos 
with opposite hypercharge. 



3.2 Grand Unification 

Before writing the Lagrangian in detail, I would like to point out that there is an interesting 
conclusion that follows from the quantum number assignments of the new particles that we have 
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Figure 5: The anomaly cancellation that requires two doublets of Higgs fields in the MSSM. 



introduced to make the Standard Model supersymmetric. 



An attractive feature of the Standard Model is that the quarks and leptons of each generation fill 
out multiplets of the simple gauge group SU{5). This suggests a very beautiful picture, called grand 
unification, in which SU{5), or a group such as S'O(IO) or Eq for which this is a subgroup, is the 
fundamental gauge symmetry at very short distances. This unified symmetry will be spontaneously 
broken to the Standard Model gauge group SU{3) x SU{2) x C/(l). 

For definiteness, I will examine the model in which the grand unified symmetry group is SU{5). 
The generators of SU(5) can be represented as 5 x 5 Hermitian matrices acting on the 5-dimensional 
vectors in the fundamental representation. To see how the Standard Model is embedded in SU{5), 
it is convenient to write these matrices as blocks with 3 and 2 rows and columns. Then the Standard 
Model generators can be identified as 



SU{3) 







SU{2) 







ctV2 




(123) 



In these expressions, is an SU{3) generator, cr°/2 is an SU{2) generator, and all of these matrices 
are normalized to tr[T'^T^] = l'^'^^- We should identify the last of these matrices with ■y/3/5 Y. 

The symmetry-breaking can be caused by the vacuum expectation value of a Higgs field in the 
adjoint representation of SU{5). The expectation value 



^1 



(124) 



commutes with the generators in (123) and fails to commute with the off-diagonal generators. So 
this vacuum expectation value gives mass to the off-diagonal generators and breaks the gauge group 
to SU{3) X SU{2) X [/(I). 

Matter fermions can be organized as left-handed Weyl fermions in the SU{5) representations 5 
and 10. The 5 is the conjugate of the fundamental representation of SU{5); the 10 is the antisym- 
metric matrix with two 5 indices. 
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(125) 



It is straightforward to check that each entry listed has the quantum numbers assigned to that 
field in the Standard Model. To compute the hypercharges, we act on the 5 with (—1) times the 



hypercharge generator in ( 123 1, and we act on the 10 with the hypercharge generator on each index. 
This gives the standard results, for example, Y = +| for the d and ^=^5 + 5— g 



for u and d. 



22 



The SU (5) covariant derivative is 

Vm = (dm - iguAiT^) , (126) 

where gu is the SU{5) gauge coupling. There is only room for one value here. So this model predicts 
that the three Standard Model gauge couplings are related by 

33 = 92 = .91 = gu , (127) 

where 

53 = ffs 52 = g 91 = \J~^9' ■ (128) 

Clearly, this prediction is not correct for the gauge couplings that we measure in particle physics. 

However, there is a way to save this prediction. In quantum field theory, coupling constants 
are functions of length scale and change their values significantly from one scale to another by 
renormalization group evolution. It is possible that the values of g' , g, and gs that we measure could 



evolve at very short distances into values that obey (127). 

I will now collect the formulae that we need to analyze this question. Let 

a. = I (129) 

for i = 1,2,3. The one-loop renormalization group equations for gauge couplings are 

dgt h 3 dai b, ^ 

dlogQ^-(4^^^ dlogQ^~Wf'- ^ ^ 

For C/(l), the coefficient bi is 

'eI^z-^e^^ (131) 



^1 3^5'^ 3^5' 

/ b 



where the two sums run over the multiplets of left-handed Weyl fermions and complex- valued bosons. 



The factors fr^ 

are the squares of the U{1) charges defined by p23|. For non-Abelian groups, the 



5 

expressions for the b coefficients are 



&--^C2(G)-|^C(r;)-^^C(r,) , (132) 



3 3^ 3 ^ 



where C2{G) and C(r) are the standard group theory coefficients. For SU{N), 

C2{G) = C{G) = N , C{N) = ^ . (133) 



The solution of the renormalization group equation (1301 is 



a~\Q)^a~\M)~^log^ . (134) 
Now consider the situation in which the three couplings gi become equal at the mass scale Mu, 



the mass scale of SU{5) symmetry breaking. Let au be the value of the ai at this scale. Using ( 134 1 
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we can then determine the Standard Model couphngs at any lower mass scale. The three ai{Q) are 
determined by two parameters. We can eliminate those parameters and obtain the relation 



(1 + B)a2^ - Bal 



where 



B = 



h - &2 



62-61 

The values of the ai are known very accurately ai Q = niz [38]: 
Qf-i = 8.50 ±0.14 = 29.57 ±0.02 a7 



59.00 ±0.02 



Inserting these values into (1351, we find 

B = 0.716 ±0.005 ±0.03 . 



(135) 
(136) 

(137) 
(138) 



In this formula, the first error is that propagated from the errors in (137) and the second is my 
estimate of the systematic error from neglecting the two- loop renormalization group coefficients and 
other higher-order corrections. 



We can compare the value of B in (1381 to the values of (136) from different models. The 



hypothesis that the three Standard Model couplings unify is acceptable only if the gauge theory 



that describes physics between mz and Mu gives a value of B consistent with (138). The minimal 
Standard Model fails this test. The values of the bi are 



62 - 



61 = 
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(139) 



where Ug is the number of generations and rih is the number of Higgs doublets. Notice that Ug cancels 



out of (136). This is to be expected. The Standard Model fermions form complete representations 
of SU{5), and so their renormalization effects cannot lead to differences among the three couplings. 
For the minimal case rih = 1 we find B = 0.53. To obtain a value consistent with (138), we need 
nil = 6. 

We can redo this calculation in the minimal supersymmetric version of the Standard Model. First 



of all, we should rewrite ( 132 ) for a supersymmetric model with one vector supermultiplet, containing 



a vector and a Weyl fermion in the adjoint representation, and a set of chiral supermultiplets indexed 



by k, each with a Weyl fermion and a complex boson. Then ( |132[ ) becomes 



h = yC2(G) - \C2{G) - 



(140) 



The formula (131) undergoes a similar rearrangement. Inserting the values of the C{rk) for the 
fields of the Standard Model, we find 

63 = 9 - 2ng 

62 = 6 - 2n„ - l^Uh 



61 = — 2n„ Uh 

1 » 10 



(141) 
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Figure 6: Prediction of the SU{3) gauge coupling from the electroweak coupling constants using grand 
unification, in the Standard Model and in the MSSM. 



For the minimal Higgs content rih — 2, this gives 



B = ^ = 0.714 (142) 



in excellent agreement with ( 138 1 



In Fig. |6] I show the unification relation pictorially. The three data points on the the left of 



the figure represent the measured values of the three couplings (1371. Starting from the values of 



ai and a2, we can integrate (1301 up to the scale at which these two couplings converge. Then we 
can integrate the equation for back down to Q = mz and see whether the result agrees with the 
measured value. The lower set of curves presents the result for the Standard Model with = 1. 
The upper set of curves shows the result for the supersymmetric extension of the Standard Model 
with n/i — 2. This choice gives excellent agreement with the measured value of a^. 

Actually, I slightly overstate the case for supersymmetry by ignoring two-loop terms in the 
renormalization group equations, and also by integrating these equations all the way down to 
even though, from searches at high-energy colliders, most of the squarks and gluinos must be heavier 
than 300 GeV. A more accurate prediction of as(mz) from the electroweak coupling constants gives 
a slightly higher value, 0.13 instead of 0.12. However, these corrections could easily be compensated 
by similar corrections to the upper limit of the integration, following the details of the particle 
spectrum at the grand unification scale. For a more detailed formal analysis of these corrections, 
see [3^, and for a recent evaluation of their effects, see HO]. It remains a remarkable fact that the 
minimal supersymmetric extension of the Standard Model is approximately compatible with grand 
unification 'out of the box', with no need for further model-building. 
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3.3 Construction of the Lagrangian 



Now I would like to write the full Lagrangian of the minimal supersymmetric extension of the 
Standard Model, which I will henceforth call the MSSM. 

The kinetic terms and gauge couplings of the MSSM Lagrangian are completely determined by 
supersymmetry, the choice of the gauge group SU (3) x SU{2) x U{1), and the choice of the quantum 



numbers of the matter fields. The Lagrangian is a sum of terms of the forms (82 1 and (84 1. Up to 



this point, the only parameters that need to be introduced are the gauge couplings gi, g2, and 173. 

Next, we need a superpotential W . The superpotential is the source of nonlinear fermi on-scalar 
interactions, so we should include the appropriate terms to generate the Higgs Yukawa couplings 
needed to give mass to the quarks and leptons. The appropriate choice is 

Wy = vid'Hdc^e^pQ^p + yl^e^Hdc^e^pL^p - fJu'E^c^^c^pQ^p . (143) 



The notation for the quark and lepton multiplets is that in (1201; the indices i,j = 1,2,3 run over 
the three generations. The indices a, (3 = 1,2 run over SU(2) isospin indices. Notice that the first 
two terms require a Higgs field Hd with Y = while the third term requires a Higgs field iJ„ with 
y = |. If we leave out one of the Higgs multiplets, some quarks or leptons will be left massless. 
This is the second argument that requires two Higgs fields in the MSSM. 



I have written ( 143 ) including the most general mixing between left- and right-handed quarks 
and leptons of different generations. However, as in the minimal Standard Model, we can remove 
most of this flavor mixing by appropriate fleld redefinitions. The coupling constants yd, Ue, Uu are 
general 3x3 complex-valued matrices. Any such matrix can be diagonalized using two unitary 
transformations. Thus, we can write 

yd = WdYdV} ye = WeYeV^ yu = WuYuVl , (144) 

with Wq and Vq 3 x 3 unitary matrices and Yq real, positive, and diagonal. The unitary transfor- 
mations cancel out of the kinetic energy terms and gauge couplings in the Lagrangian, except that 
the W boson coupling to quarks is transformed 



gu'u 



ta"dM^+ ^ gu^a^{VlYd)dWX . (145) 



From this equation, we can identify (V^Vu) = Vckm, the Cabibbo-Kobayashi-Maskawa weak inter- 
action mixing matrix. The Lagrangian term ( |143[ ) thus introduces the remaining parameters of the 
Standard Model, the 9 quark and lepton masses (ignoring neutrino masses) and the 4 CKM mixing 



angles. The field redefinition (144) can also induce or shift a QCD theta parameter, so the MSSM, 
like the Standard Model, has a strong CP problem that requires an axion or another model-building 
solution |1T] . 

There are several other terms that can be added to W. One possible contribution is a pure Higgs 
term 

Wf, = -iJ-Hdaea^Hufi . (146) 

The parameter ^ has the dimensions of mass, and consequently this mu term provides a supersym- 
metric contribution to the masses of the Higgs bosons. Because this term is in the superpotential, it 
does not receive additive raditive corrections. Even in a theory that includes grand unification and 
energies scale of the order of 10^^ GeV, we can set the parameter /Lt to be of order 100 GeV without 
finding this choice affected by large quantum corrections. We will see in Section 4.2 that the mu 
term is needed for phenomenological reasons. If /i = 0, a Higgsino state will be massless and should 
have been detected already in experiments. It is odd that a theory whose fundamental mass scale 
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is the grand unification scale should require a parameter containing a weak interaction mass scale. 
I will present some models for the origin of this term in Section 3.5. 

At this point, we have introduced two new parameters beyond those in the Standard Model. One 
is the value of fi. The other is the result of the fact that we have two Higgs doublets in the model. 
The ratio of the Higgs vacuum expectation values 

{Hu) I {Hd)=i&np (147) 

will appear in many of the detailed predictions of the MSSM. 

There are still more superpotential terms that are consistent with the Standard Model gauge 
symmetry and quantum numbers. These are 

Wfi = rjieijkUidjdk + md^apLaQp 

+ r]4eai3LaHup ■ (148) 

Here i,j,k are color indices, a, (3 are isospin indices, and arbitrary generation mixing is also possible. 
These terms violate baryon and lepton number through operators with dimensionless coefficients. 
In constructing supersymmetric models, it is necessary either to forbid these terms by imposing 
appropriate discrete symmetries or to arrange by hand that some of the dangerous couplings are 
extremely small 

If baryon number B and lepton number L are conserved in a supersymmetric model, this model 
respects a discrete symmetry called R-parity, 

i? = (-1)3^+^+2^ . (149) 

Here (SB) is quark number and J is the spin of the particle. This quantity is constructed so that 
i? = +1 on the particles of the Standard Model (including the Higgs bosons) and i? = — 1 on their 
supersymmetry partners. R acts differently on particles of different spin in the same supermultiplet, 
so R-parity is a discrete subgroup of a continuous R-symmetry. 

In a model with grand unification, there will be baryon number and lepton number violation, 
and so B and L cannot be used as fundamental symmetries. However, we can easily forbid most 



of the superpotential terms ( 148 1 by introducing a discrete symmetry that distinguishes the field 
Hd from the lepton doublets Li. A similar strategy can be used to forbid the first, 3-quark, term. 
With these additional discrete symmetries, the MSSM, including all other terms considered up to 
this point, will conserve R-parity. 



3.4 The Lightest Supersymmetric Particle 

If R-parity is conserved, the lightest supersymmetric particle will be absolutely stable. This 
conclusion has an important implication for the relation of supersymmetry to cosmology. If a 
supersymmetric particle is stable for a time longer than the age of the universe, and if this particle 
is electrically neutral, that particle is a good candidate for the cosmic dark matter. In Sections 6.3 
and 6.4, I will discuss in some detail the properties of models in which the lightest Standard Model 
superpartner is the dark matter particle. 

However, this is not the only possibility. Over times much longer than those of particle physics 
experiments — minutes, years, or billions of years — we need to consider the possibility that the lightest 
Standard Model superpartner will decay to a particle with only couplngs of gravitational strength. 
Complete supersymmetric models of Nature must include a superpartner of the graviton, a spin-| 
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particle called the gravitino. In a model with exact supersymmetry, the gravitino will be massless, 
but in a model with spontaneously broken supersymmetry, the gravitino acquires a mass through 
an analogue of the Higgs mechanism. If the supersymmetry breaking is induced by one dominant 
i^-term, the value of this mass is [33^ 

Stt (F) , , 

^3/2 = ^—. (150) 

This expression is of the same order of magnitude as the expressions for Standard Model superpartner 
masses that I will give in Section 3.6. In string theory and other unified models, there may be 
additional Standard Model singlet fields with couplings of gravitation strength, called moduli, that 
might also be light enough that long-lived Standard Model superpartners could decay to them. 

Supersymmetric models with R-parity conservation and dark matter, then, divide into two 
classes, according to the identity of the lightest supersymmetric particle — the LSP. On one hand, 
the LSP could be a Standard Model superpartner. Cosmology requires that this particle is neutral. 
Several candidates are available, including the fermionic partners of the photon, Z^, and neutral 
Higgs bosons and the scalar partner of one of the neutrinos. In all cases, these particles will be 
weakly interacting; when they are produced at high-energy colliders, they should not make signals 
in a particle detector. On the other hand, the LSP could be the gravitino or another particle with 
only gravitational couplings. In that case, the lightest Standard Model superpartner could be a 
charged particle. Whether this particle is visible or neutral and weakly interacting, its decay should 
be included in the phenomenology of the model. 



3.5 Models of Supersymmetry Breaking 

There is still one important effect that is missing in our construction of the MSSM. The terms that 
we have written so far preserve exact supersymmetry. A fully supersymmetric model would contain a 
massless fermionic partner of the photon and a charged scalar particle with the mass of the electron. 
These particles manifestly do not exist. So if we wish to build a model of Nature with supersymmetry 
as a fundamental symmetry, we need to arrange that supersymmetry is spontaneously broken. 

From the example of spontaneous symmetry breaking in the Standard Model, we would expect to 
do this by including in the MSSM a field whose vacuum expectation value leads to supersymmetry 
breaking. This is not as easy as it might seem. To explain why, I will first present some models of 
supersymmetry breaking. 

The simplest model of supersymmetry breaking is the O'Raifeartaigh model [44] , with three 
chiral supermultiplets (f)o, 4>i, 4>2 interacting through the superpotential 

W^X(l)a + m0i02 + g<f>o<l>l ■ (151) 

This superpotential implies the F = conditions 

= F*^X + gq^l 

= F* = TO02 + 2g(j)o(j)i 

Q = F; ^ m<t>i (152) 

The first and third equations contradict one another. It is impossible to satisfy both conditions, 
and so there is no supersymmetric vacuum state. This fulfils the condition for spontaneous super- 
symmetry breaking that I presented in Section 2.2. 
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This mechanism of supcrsymmetry breaking has an unwanted corollary. Because one combina- 



tion of the scalar fields appears in two different constraints in (1521, there must be an orthogonal 
combination that does not appear at all. This means that the F-term potential Vp has a surface of 
degenerate vacuum states. To see this explicitly, pick a particular vacuum solution 

00 - (/-i - 02 = . (153) 

and expand the potential Vp about this point. There are 6 real- valued boson fields with masses 



, , m , m , \/m2 - 2Xg , y^m^ + 2\g . (154) 
These six fields do not pair into complex-valued fields; that is already an indication that supcrsym- 



metry is broken. The fermion mass term in (38) gives one Dirac fermion mass m and leaves one 



Weyl fermion massless. This massless fermion is the Goldstone particle associate with spontaneous 
supcrsymmetry breaking. 

A property of these masses is that the sum rule for fermion and boson masses 

str[m2] ^^m) ~^ml ^ Q (155) 

remains valid even when supcrsymmetry is broken. This sum rule is the coefficient of the one-loop 
quadratic divergence in the vacuum energy. Since supcrsymmetry breaking does not affect the ul- 
traviolet structure of the theory, this coefficient must cancel even if supcrsymmetry is spontaneously 
broken |45j . In fact, if Q is a conserved charge in the model, the sum rule is valid in each charge 
sector Q = q: 

strjm^] = . (156) 

In the O'Raifeartaigh model, supcrsymmetry is spontaneously broken by a nonzero expectation 
value of an F term. It is also possible to break supcrsymmetry with a nonzero expectation value of 
a D term. The I?-term potential Vd typically has zeros. For example, in an SU{S) supersymmetric 
Yang-Mills theory, 

^d^I{t. '^^^"•^ - E ^^"^^ j (157) 

and it is easy to find solutions in which the terms in parentheses sum to zero. However, it is not 
difficult to arrange a Vp such that the solutions of the F — conditions do not coincide with the 
solutions of the D = conditions. This leads to spontaneous symmetry breaking, again with the 



sum rule ( 156 1 valid at tree level. 



Unfortunately, the sum rule (1561 is a disaster for the prospect of finding a simple model of 
spontaneously broken supcrsymmetry that extends the Standard Model. For the charge sector of 
the d squarks, we would need all down-type squarks to have masses less than 5 GeV. For the charge 
sector of the charged leptons, we would need all sleptons to have masses less than 2 GeV. 



3.6 Soft Super symmetry Breaking 

The solution to this problem is to construct models of spontaneously broken supcrsymmetry using 
a different strategy from the one that we use for electroweak symmetry breaking in the Standard 
Model. To break electroweak symmetry, we introduce a Higgs sector whose mass scale is the same 
as the scale of the fermion and gauge boson masses induced by the symmetry breaking. To break 
supcrsymmetry, however, we could introduce a new sector at a much higher mass scale, relying on a 
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weak coupling of the new sector to the Standard Model particles to communicate the supersymmetry 
breaking terms. In principle, a weak gauge interaction could supply this coupling. However, the 
default connection is through gravity. Gravity and supergravity couple to all fields. It can be shown 
that supersymmetry breaking anywhere in Nature is communicated to all other sectors through 
supergravity couplings. 

We are thus led to the following picture, which produces a phenomenologically reasonable su- 
persymmetric extension of the Standard Model: We extend the Standard Model fields to supersym- 
metry multiplets in the manner described in Section 3.1. We also introduce a hidden sector with 
no direct coupling to quark, leptons, and Standard Model gauge bosons. Supersymmetry is spon- 
taneously broken in this hidden sector. A weak interaction coupling the two sectors then induces a 
supersymmetry-breaking effective interaction for the Standard Model particles and their superpart- 
ners. If A is the mass scale of the hidden sector, the supersymmetry breaking mass terms induced 
for the Standard Model sector are of the order of 

— ^-tJ: (158) 

MM ^ ^ 

where M is the mass of the particle responsible for the weak connection between the two sectors. 
M is called the messenger scale. By default, the messenger is supergravity. Then M = mpi and 
A ~ 10^^ GeV. In this scenario, the superpartners acquire masses of the order of the parameter m 



in (1581 



It remains true that the quarks, leptons, and gauge bosons cannot obtain mass until SU{2) x 
U{1) is broken. It is attractive to think that the symmetry-breaking terms that give mass to the 
superpartners cause SU{2) x C/(l) to be spontaneously broken, at more or less the same scale. I 
will discuss a mechanism by which this can happen in Section 6.1. The weak interaction scale would 
then not be a fundamental scale in Nature, but rather one that arises dynamically from the hidden 
sector and its couplings. 

The effective interaction that are generated by messenger exchange generally involve simple 
operators of low mass dimensions, to require the minimal number of powers of M in the denominator. 
These operators are soft perturbations of the theory, and so we say that the MSSM is completed by 
including soft supersymmetry-breaking interactions. 

However, the supersymmetry-breaking terms induced in this model will not include all possible 
low-dimension operators. Since these interactions arise by coupling into a supersymmetry theory, 
they are formed by starting with a supersymmetric effective action and turning on F and D ex- 
pectation values as spurious. Only a subset of the possible supersymmetry-breaking terms can be 
formed in this way |46j . By replacing a superfield <i> by 6^c9 (F), we can convert 

J d^0K{<^>,<j)) m20t0 
J d^e f{^)W'^cW ^ mX^cX 

d^W{<P,<j)) B(j)^ + A(f>^ (159) 



However, as long as the (j) theory is renormalizable, we cannot generate the terms 

mV'^cV' , C(j)*4>^ , (160) 

by turning on expectation values for F and D fields. Thus, we cannot generate supersymmetry-break- 
ing interactions that are mass terms for the fermion field of a chiral multiplet or non-holomorphic 
cubic terms for the scalar fields. 
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There is another difficulty with terms of the form ( 160 1 . In models with Standard Model singlet 
scalar fields, which typically occur in concrete models, these two interactions can generate new 
quadratic divergences when they appear in loop diagrams [46J. 

Here is the most general supersymmetry-breaking effective Lagrangian that can be constructed 
following the rule just given that is consistent with the gauge symmetries of the Standard Model: 



ft 



Ta 



-Mf\ff--m,X, 

-{AdUddHda^apQp 



cXf 



AeUeeHdaf-aflLp 
AuyuUHua<^apQl3 - B ^Hda^alsHufi) - h.C. 



(161) 



I have made the convention of scaling the A terms with the corresponding Yukawa couplings and 
scaling the B terms with fj,. The parameters A and B then have the dimensions of mass and are 
expected to be of the order of m in ( 158 1 . 

For most of the rest of these lectures, I will represent the effects of the hidden sector and 
supersymmetry breaking simply by adding (1611 to the supersymmetric Standard Model. I will 
then consider the MSSM to be defined by 



jw ■ 



^soft 



(162) 



combining the pieces from (82l, (84), (1431, (1461, and (161|. 



There are two problems with this story. The first is the /i term in the MSSM superpotential. 
This a supersymmetric term, and so fj. can be arbitrarily large. To build a successful phenomenology 
of the MSSM, however, we need to have n of the order of the weak scale. Ideally, fi should be 
parametrically equal to (1581. 

There are simple mechanisms that can solve this problem. A fundamental theory that leads 
to the renormalizable Standard Model at low energies can also contain higher-dimension operators 
suppressed by the high-energy mass scale. Associate this scale with the messenger scale. Then a 
supersymmetric higher-dimension operator in the superpotential 



(163) 



leads to a /X term if S acquires a vacuum expectation value. If 5 is a hidden sector field, we could 
find gZ] 



^ M M ' 

A supersymmetric higher dimension contribution to the Kahler potential 

(fe^<i>^HdH^ 



(165) 



leads to a /X term if $ acquires a vacuum expectation value in its F term. If <I> is a hidden sector 
field, we could find |15] 

In models with weak-coupling dynamics, higher-dimension operators are associated with the string or 
Planck scale; then, these mechanisms work most naturally if supergravity is the mediator. However, 
it is also possible to apply these strategies in models with strong-coupling dynamics in the hidden 
sector at an intermediate scale. 
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Figure 7: A dangerous contribution to K-K mixing involving gluino exchange and flavor mixing in the 
squark mass matrix. 



Generating the /i term typically requires breaking all continuous R-symmetries of the model. 
This is unfortunate, because an R-symmetry might be helpful phenomenologically, for example, to 
keep gaugino masses small while allowing sfermion masses to become large, or because it might be 
difficult to break an R-symmetry using a particular explicit mechanism of supersymmetry breaking. 
In this case, it is necessary to add Standard Model singlet fields to the MSSM to allow all gaugino 
and Higgsino fields to acquire nonzero masses. Models of this type are presented in [ 49150] . 

The second problem involves the flavor structure of the soft supersymmetry breaking terms. 



In writing (161 1, I did not write flavor indices. In principle, these terms could have flavor- mixing 



that is arbitrary in structure and different from that in ( 143 ) . Then the flavor- mixing would not 



be transformed away when ( 143 1 is put into canonical form. However, flavor-mixing from the soft 
supersymmetry breaking terms is highly constrained by experiment. Contributions such as the one 
shown in Fig. [t] give contributions to K'^ , , and mixing, and to r — > /i7 and /i 67, that 
can be large compared to the measured values or limits. Theories of the origin of the soft terms 
in models of supersymmetry breaking should address this problem. For example, the models of 
gauge-mediated [52 and anomaly-mediated [53154] supersymmetry breaking induce soft terms that 
depend only on the SU{2) x U{1) quantum number and are therefore automatically diagonal in 
flavor. A quite different solution, based on a extension of the MSSM with a continuous R-symmetry, 
is presented in |5IJ. 

If I assume that the soft supersymmetry-breaking Lagrangian is diagonal in flavor but is otherwise 
arbitrary, it introduces 22 new parameters. With arbitrary flavor and CP violation, it introduces 
over 100 new parameters. This seems a large amount of parameter freedom. I feel that it is not 
correct, though, to think of these as new fundamental parameters in physics. The soft Lagrangian 
is computed from the physics of the hidden sector, and so we might expect that these parameters 
are related to one another as a part of a theory of supersymmetry breaking. Indeed, the values of 
these parameters are the essential data from which we will infer the properties of the hidden sector 
and its new high energy interactions. 

If supersymmetry is discovered at the weak interaction scale, it will be a key problem to measure 
the coefficients in the soft Lagrangian and to understand their pattern and implications. Most of 
my discussion in the next two sections will be devoted to the question of how the soft parameters 
can be determined from data at the LHC and ILC. 
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4 The Mass Spectrum of the MSSM 



4.1 Sfermion Masses 



Our first task in this program is to ask how the parameters of the MSSM Lagrangian are reflected 
in the mass spectrum of the superparticles. The relation between the MSSM parameters and the 
particle masses is surprisingly complicated, even at the tree level. For each particle, we will need 
to collect all of the pieces of the Lagrangian ( 162 1 that can contribute to the mass term. Some of 
these will be direct mass contributions; others will contain Higgs fields and contribute to the masses 
when these fields obtain their vacuum expectation values. In this discussion, and in the remainder 
of these lectures, I will ignore all flavor-mixing. 

Begin with the squark and slepton masses. For light quarks and leptons, we can ignore the 
fermion masses and Higgs couplings. Even with this simplification, though, there are two sources 
for the scalar masses. One is the soft mass term 

Lsoft = -M/l/P . (167) 

The other comes from the Z?-term potential. The SU{2) and C/(l) potentials contain the cross terms 
between the Higgs field and sfermion field contributions 

+ ^ • 2 • {-IhIh, + \HiH^) ■ irrf) . (168) 

To evalute this expression, we must insert the vacuum expectation values of the two Higgs fields. 
In terms of the angle (3 defined in (147), these are 

w,>-(^.i,) (i/.)^(*"r''), (169) 

where v = 246 GeV so that rriw = gv/2. 



Inserting the Higgs vevs into the potential (1681, we find 

Vd = r [^(cos^ (3 - sin' P)ig'l' - a''Y)\ f 
= r [^l±p^ cos 2/3 (/3 - slil' + Y))]f 



= r[m|cos2/3(/3-4Q)]/ . (170) 



Then, if we define 



Af = {l'~siQ)cos2l3m'z, (171) 
the mass of a first- or second-generation sfermion takes the form 

mj = + Af (172) 

when contributions proportional to fermion masses can be neglected. The D-term contribution can 
have interesting effects. For example, SU{2) invariance of M| implies that 

m2(e) - m'(j^) = |cos2/3| to| > . (173) 
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For some choices of parameters, the measurement of this mass difference is a good way to determine 
tan/3 ISg. 

For third-generation fermions, the contributions to the mass term from Yukawa coupHngs and 

from A terms can be important. For the b and b, these contributions come from the terms in the 
effective Lagrangian 

\F,\' + \F^\' - \y, {H',)b\' + \yJb{Hi) p = mlm^ + m 

l-FffdP = {- 11 {H^d))* {ybbb) + h.c. = - fimb tan f3 bb + h.c. 
-Lsoft = Abyb (H^) lb = Atmblb . (174) 
In all, we find a mass matrix with mixing between the two scalar partners of the b quark, 

T 



{b* f)Ml(^^ 



(175) 



with 

+ Afc + ml mbiAb- tan /3) 



{^biAb- fitan/3) x2^A-^-2 ) (176) 



The mass matrix for r, t has the same structure. For t, t, replace tan /3 by cot /?. 

The mixing terms in the mass matrices of the third-generation sfermions often play an important 
role in the qualitative physics of the whole SUSY model. Because of the mixing, one sfermion 
eigenstate is pushed down in mass. This state is often the lightest squark or even the lightest 
superparticle in the theory. 

4.2 Gaugino and Higgsino Masses 

In a similar way, we can compute the mass terms for the gauginos and Higgsinos. Since the 
gauginos and Higgsino have the same quantum numbers after SU{2) x U{1) breaking, they will mix. 
We have seen in Section 2.4 that this mixing plays an essential role in the working of the Higgs 
mechanism in the limit where soft supersymmetry breaking terms are turned off. 

The charged gauginos and Higgsinos receive mass from three sources. First, there is a soft SUSY 
breaking term 

- Lso/t = m2W^'^cw^ . (177) 
The /i superpotential term contributes 

-Lw^ ^l'Ki^cht . (178) 

The gauge kinetic terms contribute 

- L = ^j=^m) ^'^^ch-, + {Hi) w-^cht) (179) 



Inserting the Higgs field vevs from (169), we find the mass term 



with 



{w-^ h^^)cmci~+], (180) 



, 7712 V^TOvySm/? 1 , . 

' Vzniw COS p ^ ' 
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The mass matrix for neutral gauginos and Higgsinos also receives contributions from these three 
sources. In this case, all four of the states 



(182) 



have the same quantum numbers after SU(2) x U{1) breaking and can mix together. The mass 
matrix is 

/mi -mzCfjSw mzsps^^ 

TO2 mzCfsCw -mzspcw 

-mzCfjSw mzCfsCw 

V mzSfjSw -mzspCw -fJ- 



rriN 



(183) 



The mas s eige nstates in these systems are referred to collectively as charginos and neutralinos. 
The matrix (1831 is complex symmetric, so it can be diagonalized by a unitary matrix Vo[^ 



TUN = V^DnV^ 



(184) 



I will denote the neutralinos as Nf , i = 1, . . . , 4, in order of mass with the lightest. Elsewhere 
in the literature, you will see these states called or The mass eigenstates are related to the 
weak eigenstates by the transformation 



Vo 



N2 



(185) 



Note that the diagonal matrix in (184) may have negative or complex-valued elements. If that 



is true, the physical fermion masses of the Ni are the absolute values of the corresponding elements 
of Dat. The phases will appear in the three-point couplings of the Ni and can lead to observable 
interference effects. Complex phases in would provide a new source of CP violation. 



The chargino mass matrix ( 181 1 is not symmetric, so in general it is diagonalized by two unitary 
matrices 



mc = VlDcV_ 



(186) 



I will denote the charginos as 



1,2, in order of mass with the lighter. Elsewhere in the 



Hterature, you will see these states called xf or . The mass eigenstates are related to the weak 
eigenstates by the transformation 



h+ 



ct 



5r 



(187) 



It should be noted that /i are must be nonzero. If /i = 0, the determinant of (1831 vanishes and 



so the lightest neutralino must be massless. This neutralino will also have a large Higgsino content 
and thus an order- 1 coupling to the . It is excluded by searches for an excess of invisible 
decays and for Z° NiN2- The condition /i = also implies that the lightest chargino has a mass 
below the current limit of about 100 GeV. 

Often, one studies models for which mi, m2, and /i are all large compared to mw and mz- The 
off-diagonal elements that mix the gaugino and Higgsino states are of the order of mw and mz- 



■l-Note that this formula is different from that which diagonalizes a Hermitian matrix. A detailed discussion of the 
diagonalization of mass matrices appearing in SUSY can be found in the Appendix of [56) . 
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Figure 8: Masses of the four neutralinos and two charginos along a line in the SUSY parameter space on 
which m{Ni) — 100 GeV while the parameter fi moves from large negative to large positive values. The 
parameter mi is set to rrii — 0.5m2. Note the approximate degeneracies in the extreme limits of the gaugino 
and Higgsino regions. 



Thus, if the scale of masses generated by the SUSY breaking terms is large, the mixing is small and 
the individual eigenstates are mainly gaugino or mainly Higgsino. However, there are two distinct 
cases. The first is the gaugino region, where mi,r7i2 < IfJ-l- In this region of parameter space, the 
lightest states A^i, Ci arc mainly gaugino, while the heavy neutralinos and charginos are mainly 
Higgsino. In the Higgsino region, mi,m2 > the situation is reversed and Ni, Ci are mainly 
Higgsino. In this case, the two lightest neutralinos are almost degenerate. In Fig. [8] I show the 
mass eigenvalues as a function of the mass matrix parameters along a line in the parameter space 
on which the Ni has a fixed mass of 100 GeV. As we will see in Section 6.4, the exact makeup of 
the lightest ncutralino as a mixture of gaugino and Higgsino components is important to the study 
of supersymmetric dark matter. 

To summarize this discussion, I present in Fig. [9] the complete spectrum of new particles in the 
MSSM at a representative point in its parameter space. Notice that the third-generation sfermions 
are split off from the others in each group. Note also that the parameter point chosen is in the 
gaugino region. The lightest superparticle is the A^i. I will discuss the spectrum of Higgs bosons in 
Section 6.2. 



4.3 Renormalization Group Evolution of MSSM Parameters 

The spectrum shown in Fig. [9] appears to have been generated by assigning random values to 
the soft SUSY breaking parameters. But, actually, I generated this spectrum by making very simple 
assumptions about the relationships of the soft parameters, at a high energy scale. Specifically, I 
assumed that the soft SUSY breaking gaugino masses and (separately) the sfermion masses were 
equal at the scale of grand unification. The structure that you see in the figure is generated by the 
renormalization group evolution of these parameters from the grand unification scale to the weak 
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Figure 9: Illustrative spectrum of supersymmetric particles. The columns contain, from the left, the Higgs 
bosons, the four neutralinos, the two charginos, the charged sleptons, the sneutrinos, the down squarks, and 
the up squarks. The gluino, not shown, is at about 800 GeV. 



scale. 

The renomalization group (RG) evolution of soft parameters is likely to play a very important 
role in the interpretation of measurements of the SUSY particle masses. Essentially, after measuring 
these masses, it will be necessary to decode the results by running the effective mass parameters up 
to a higher energy at which their symmetries might become more apparent. The situation is very 
similar to that of the Standard Model coupling constants, where a renormalization group analysis 



told us that the apparently random values ( 137 1 for the coupling constants at the weak scale actually 



corresponds to a unification of couplings at a much higher scale. 

In this section, I will write the most basic RG equations for the soft gaugino and sfermion masses. 
One further effect, which involves the Yukawa couplings and is important for the third generation, 
will be discussed later in Section 6.1. 

The RG equation for the gaugino masses is especially simple. This is because both the gaugino 



masses and the gauge couplings arise from the superpotential term (115), with the supersymmetry 



breaking terms arising as shown in (1591. As I have already noted, this F-term receives a radiative 
correction proportional to the /3 function as a consequence of the trace anomaly |34l36j . The correc- 
tions are the same for the gauge boson field strength and the gaugino mass. Thus, if gaugino masses 
and couplings are generated at the scale M, they have the relation after RG running to the scale Q: 

mi{Q) _ ai{Q) 



TOi(M) a,{M) 



(188) 



If the F term that generates the soft gaugino masses is an SU{5) singlet, the soft gaugino masses 
will be grand-unified at M. Then, running down to the weak scale, they will have the relation 

mi : 1712 : — ai : a2 ■ = 0.5 : 1 : 3.5 . (189) 
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Figure 10: Diagrams that generate the soft mass parameters in gauge mediated supersymmetry breaking 
(a.) gaugino masses; (b.) sfermion masses. 
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Figure 11: Diagrams that generate the renormalization group evolution of the soft sfermion mass parameters 
Mj. 



This relation of soft gaugino masses is known as gaugino unification. 

There are other models of the soft gaugino masses that also lead to gaugino unification. In 
gauge-mediated SUSY breaking, the dynamics responsible for SUSY breaking occurs at a scale much 
lower than the scale associated with mediation by supergravity. At this lower scale Mg (for example, 
1000 TeV), some heavy particles with nontrivial SU{3) x SU{2) x U{1) quantum numbers acquire 
masses from SUSY breaking. These fields then couple to gauginos and generate SUSY breaking 
masses for those particles through the diagram shown in Fig. [Toj[a). The heavy particles must fall 
into complete SU{5) representations; otherwise, the coupling constant renormalization due to these 
paxticles between Mg and the grand unification scale would spoil the grand unification of the gauge 



couplings. Then the diagram in Fig. 10 'a) generates soft gaugino masses proportional to a{Mg). 



Running these parameters down to the weak scale, we derive the relation (1891 from this rather 
different mechanism. 

Now let us turn to the RG running of soft scalar masses. In principle, there are two contributions, 
one from the RG rescaling of the soft mass term M| and one from RG evolution generating M| 
from the gaugino mass. The Feynman diagrams that contribute to the RG coefficients are shown in 
Fig. [U] The two one-loop diagrams proportional to Mj cancel. The third diagram, involving the 
gaugino mass, gives the RG equation 



dMf _ 2 



dlogQ 



^«,(Q)C2(r,)m2(Q) , (190) 



with i — 1,2,3 and C2{ri) the squared charge in the fermion representation under the gauge 
group i. This equation leads to a positive contribution to Mj as one runs the RG evolution from 
the messenger scale down to the weak scale. The effect is largest for squarks, for which the SUSY 
breaking mass is induced from the gluino mass. 

As an example of this mechanism of mass generation, assume gaugino unification and assume 
that Mj = for all sfermions at the grand unification scale. Then the weak scale sfermion masses 
will be in the ratio 

M(i) : M(e) : M(d) : M{u) : M{d,u) : 

= 0.5 : 0.9 : 3.09 : 3.10 : 3.24 : 1 (191) 
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Figure 12: Evolution of squark and slepton masses from the messenger scale down to the weak scale, for four 
different models of supersymmetry breaking: (a.) universal sfermion masses at the grand unification scale 
Mu\ (b.) sfermion masses at Mu that depend on the SU{^) representation; (c.) universal sfermion masses 
at an intermediate scale; (d.) gauge mediation from a sector of mass about 1000 TeV. 



This model of fcrmion mass generation is called no-scale SUSY breaking. It has the danger that 
the lightest stau mass eigenstate could be lighter than than the Ni, leading to problems for dark 
matter. This problem can be avoided by RG running above the GUT scale [57]. Alternatively, it 
might actually be that the lightest Standard Model superpartner is a long-lived stau that eventually 
decays to a tau and a gravitino |58l59j . 

In gauge-mediated SUSY breaking, the diagram shown in Fig. [Toj^b) leads to the qualitatively 
similar but distinguishable formula 

a2 



M} = 2Y,al{M)C2{n)-['^] . (192) 



Each model of SUSY breaking leads to its own set of relations among the various soft SUSY 
breaking parameters. In general, the relations are predicted for the parameters defined at the mes- 
senger scale and must be evolved to the weak scale by RG running to be compared with experiment. 



Fig. 12 shows four different sets of high-scale boundary conditions for the RG evolution, and the 
corresponding evolution to the weak scale. If we can measure the weak-scale values, we could try to 
undo the evolution and recognize the pattern. This will be a very interesting study for the era in 
which superparticles are observed at high energy colliders. 
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There are some features common to these spectra that are general features of the RG evolution 
of soft parameters: 

1. The pairs of sleptons e and e can easily acquire a significant mass difference from RG evolution, 
and they might also have a different initial condition. It is important to measure the mass 
ratio m{e)/m{e) as a diagnostic of the scheme of SUSY breaking. 

2. Gaugino unification is a quantitative prediction of certain schemes of SUSY breaking. It is im- 
portant to find out whether this relation is correct or not for the real spectrum of superparticles 
in Nature. 

3. When the RG effects on the squark masses dominate the values of M| from the initial condition, 
the various species of squark have almost the same mass and are much heavier than the sleptons. 
It is important to check whether most or all squarks appear at the same threshold. 



5 The Measurement of Supersymmetry Parameters 

5.1 Measurements of the SUSY Spectrum at the ILC 

Now that we have discussed the physics that determines the form of the spectrum of superparti- 
cles, we turn to the question of how we would determine this spectrum experimentally. This is not as 
easy as it might seem. In this section, I will consider only models in which the dark matter particle 
is the A^i , and all other SUSY particles decay to the Ni . This neutral and weakly interacting particle 
would escape a collider detector unseen. Nevertheless, methods have been worked out not only to 
measure the masses of superparticles but also to determine mixing angles and other information 
needed to convert these masses to values of the underlying parameters of the MSSM Lagrangian. 

Similar methods apply to other scenarios. For example, in models in which the neutralino decays 
to a particle with gravitational interactions, one would add that decay, if it is visible, to the analyses 
that I will present. It is possible in models of this type that the lightest Standard Model superpartner 
would be a charged slepton that is stable on the time scale of particle physics experiments. That 
scenario would produce very striking and characteristic events |58j . 

Most likely, this experimental study of the SUSY spectrum will begin in the next few years 
with the LHC experiments. However, at a hadron collider like the LHC, much of the kinematic 
information on superparticle production is missing and so special tricks are needed even to measure 
the spectrum. The study of supersymmetry should be much more straightforward at an e'^e~ collider 
such as the planned International Linear Collider (ILC). For this reason, I would like to begin my 
discussion of the experiments in this section by discussing SUSY spectrum measurements at e"^e~ 
colliders. More complete reviews of SUSY measurements at linear colliders can be found in |60I61| . 

I first discuss slepton pair production, beginning with the simplest process, e+e~ — > fi^Jir and 
considering successively the production of r and e. Each step will bring in new complexities and 
will allow new measurements of the SUSY parameters. 

The process e^e~ — > Jl'^]l~ , where Jl is the partner of either the left- or right-handed /i, can 
be analyzed with the simple formulae for scalar particle-antiparticle production. The cross section 
for pair production from polarized initial electrons and positrons to final-state scalars with definite 
SU(2) X C/(l) quantum numbers is given by 

-^ = '^f3'sinH\U\ (193) 
a cos u 2s 
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Figure 13: Schematic energy distribution of final-state muons in e^e — > ^"""/i 



where 

and, in this expression, /■^ = — ij, for a, 5 = L,i?. For the initial state, a ^ L denotes the state 
e£ejj and a ^ R denotes e^ej. For the final state, b — L denotes the Jl, b = R the JI. Notice that 
this cross section depends strongly on the polarization states: 

l/abP = 1.69 e~^e\-^Ji^Ji 

= 0.42 ele+^^Ji^jT 

= 0.42 e^ej —> Tif^Ji^ 

= 1.98 e^e+ -> /7+pI~ (195) 

The angular distribution is characteristic of pair-production of a spin particle; the normalization 
of the cross sections picks out the the correct set of SU{2) x U{1) quantum numbers. 

If the smuon is light, its only kinematically allowed decay might be /i • Even if the 

smuon is heavy, if the Ni is mainly gaugino, this decay should be important. As noted above, I 
am assuming that R-parity is conserved and that the A^i is the lightest particle in the superparticle 
spectrum. Then events with this decay on both sides will appear as 

e^e~ — !■ + (missing E and p) (196) 

The spectrum of the observed muons is very simple. Since the Jl has spin 0, it decays isotropically 
in its own rest frame. In e+e~ production at a definite center of mass energy, the Jl is produced at 
a definite energy, and thus with a definite boost, in the lab. The boost of an isotropic distribution 
is a flat distribution in energy. So, the muon energy distribution should be flat, between endpoints 



determined by kinematics, as shown in the idealized Fig. 13 



The endpoint positions are simple functions of the mass of the fi and the mass of the Ni, 

(197) 

where 7 = i?cM/2m(/I), P — {1 — 4to^(/I)/_Bq,^)^/^. If we can identify both endpoint positions, we 
can solve for the two unknown masses. Figure [14] shows a simulation of the reconstructed smuon 
energy distribution from JI pair production at the ILC [63j . The high-energy edges of the distributions 
are rounded because of initial-state radiation in the e~^e~ collision. The experimenters expect to be 
able to measure this effect and correct for it. Then they should obtain values of the smuon mass to 
an accuracy of about one hundred MeV, or one part per mil. 
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Figure 14: Energy distribution of muons from e^e~ JI Ji^ at the ILC, in a simulation by Blair and 
Martyn tiiat includes realistic momentum resolution and beam effects |63| . 



A similar analysis applies to e^e 



but there are several complications. First, for the 



T system, mixing between the r and the t might be important, especially if tan/3 is large. The 
production cross sections are affected directly by the mixing. For example, to compute the pair- 
production of the lighter r mass eigenstate from a polarized initial state. 



generalize (1931 to 



da 
dcos ( 



ttq; 



where 



fui = Irr cos^ Or + Jrl sin^ Or 



we must 



(198) 



(199) 



and 6r is the mixing angle associated with the diagonalization of the r case of (1761 



Second, while the r can decay to through gauge couplings, this weak eigenstate can also 
decay to Tj^hd through terms proportional to the Yukawa coupling. Both decay amplitudes con- 
tribute to the observable decay ti tNi. With the r mixing angle fixed from the measurement of 
the cross section, the r polarization in r decays can be used to determine the mixing angles in the 



diagonalization of the neutralino mass matrix (1831 [62] 



In Fig. 15 I show the distribution of total visible energy in r — > 37r 



+ at the ILC. 

Though there is no longer a sharp feature at the kinematic endpoint, it is still possible to accurately 
determine the r mass by fitting the shape of this distribution. 

The physics of e+e~ — > e+e^ brings in further new features. In this case, there is a new 
Feyman diagram, involving t-channel neutralino exchange. The two diagrams contributing to the 
cross section for this process are shown in Fig. [16] The t-channel diagram turns out to be the more 
important one, dominating the s-channel gauge boson exchange and generating a large forward peak 



in selectron production. The cross section for e^jcj 
( fl93| , 

da 



e e is given by another generalization of 



dcos ( 



2s 



13' sin' eiJ^Riil 



where 



^RR = Jrr 



E 



(200) 



(201) 



42 



600 



400 



200 




20 40 60 80 100 
3 = prong Eg^ [GeV] 



Figure 15: Energy distribution of the three-pion system from e^e — > Tj^ at the ILC, with a r decay to 
Stt, in a simulation by Blair and Martyn that includes realistic momentum resolution and beam effects. |63j . 




Figure 16: Feynman diagrams contributing to e^e e 



with the sum running over neutralino mass eigenstates. The factor Vqu is a matrix element of the 



unitary matrix introduced in (184) 



The t-channel diagram also allows new processes such as e^ej — s- e e . Note the correlation 
of the initial-state electron and position spins with the identities of the final-state selectrons. A 
complete set of polarized cross sections for selectron pair production in e^e~ and e~e~ collisions 
can be found in [64] , 

The cross sections for chargino and neutralino pair production in e+e^ collisions are somewhat 
more complicated, but still there are interesting things to say about these processes. Chargino 
pair production is given by the Feynman diagrams shown in Fig. |17[ These diagrams are just the 
supersymmetric analogues of the diagrams for e+e^ W~^W~ . As in that process, the most 
charcteristic final states are those with a hadronic decay on one side of the event and a leptonic 
decay on the other side, for example. 



C+ i+i^N^ , Cf duN^ 



(202) 



A typical event of this kind is shown in Fig. [18] 

The chargino and neutralino production cross sections have a strong dependence on the mixing 



angles in (1841 and (186) and offer a number of strategies for the determination of these mixing 
angles. Let me present one such strategy here. Consider the reaction from a polarized initial 
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Figure 19: Contours of constant cross section for the process e^ej — > (in fb, for E^yi — 500 GeV), as 

a function of the underlying SUSY parameters |66J. The region shown is that in which the lightest chargino 
mass varies from 50 to 200 GeV. For fixed mass, the cross section increases from zero to about 150 fb 
as we move from the gaugino region into the Higgsino region. 



state e^ej — s- . Since we have an initial e^, the t-channel diagram vanishes because the right- 

handed electron does not couple to the neutrino. Now simplify the s-channel diagram by considering 
the limit of high energies, s ^ m^. In this limit, it is a good approximation to work with weak 
gauge eigenstates {B'^^W'^) rather than the mass eigenstates (7, The weak eigenstate basis 

gives a nice simplification. The initial couples only to 5°. But couple only to , so at high 
energy the s-channel diagram gets contributions only from the Higgsino components of the and 
Ci eigenstates. If we go to still higher energies, s ^ m{Cxf' , there is a further simplification. The 
cross section for h^h^ production is forward-peaked, and the cross section for hj^h"^ production is 
backward-peaked. Then, the cross section for e]^ej CfCj^ takes the form 

'^'^ ™' [|y+2l|^l + COS0)2+|F_2l|'(l-COS0)2] . (203) 



In this limit, it is clear that we can read off both of the mixing angles in (186 1 from the shape of 
this cross section. 

The use of high-energy limits simplified this analysis, but the sentivity of this cross section to 
the chargino mixing angles is not limited to high energy. Even relatively close to threshold, the 
polarized cross sections for chargino production depend strongly on the chargino mixing angles and 
can be used to determine their values. In Fig. [19] I show contours of constant cross section for 
e^ej C^Ci in the (m2,/i) plane (for tan/3 = 4 and assuming gaugino unification) The 
value of this cross section is always a good measure of whether the SUSY parameters in Nature put 
us in the gaugino or the Higgsino region of Fig. [S) 
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5.2 Observation of SUSY at the LHC 



Now we turn to supersymmetry production processes at the LHC. This subject, though more 
difficult, has immediate importance, since the LHC experiments are just about to begin. 

The reactions that produce superparticles are typically much more complicated at hadron col- 
liders than at lepton colliders. This is true for several reasons. High energy collisions of hadrons are 
intrinsically more complicated because the final states include the fragments of the initial hadrons 
that do not participate in the hard reaction. More importantly, the dominant reactions at hadron 
colliders are those that involve strongly interacting superparticles. This means that the primary 
particles are typically the heavier ones in the spectrum, which then decay in several steps. In addi- 
tion, large backgrounds from QCD obscure the signatures of supersymmetric particle production in 
many channels. 

Because of these difficulties, there is some question whether SUSY particle production can be 
observed at the LHC. However, as I will explain, the signatures of supersymmetry are still expected 
to be striking and characteristic. It is not so clear, though, to what extent it is possible to measure 
the parameters of the SUSY Lagrangian, as I have described can be done from ILC experiments. 
This is an important study that still offers much room for new ideas. 

The discovery of SUSY particles at the LHC and the measurement of SUSY parameters has been 
analyzed with simulations at a number of parameter points. Collections of interesting studies can 
be found in |63l67l68j . 

The dominant SUSY production processes at the LHC are 

gg ^99 , qq* gq-^m (204) 

These cross sections are large — tens of pb in typical cases. The values of numerous SUSY production 



cross sections at the LHC are shown in Fig. 20 |7D] 



We have seen that the squarks and gluinos are typically the heaviest particles in the supersym- 
metry spectrum. The gluinos and squarks thus will decay to lighter superparticles. Some of these 
decays are simple, e.g., 

^ qN^ . (205) 
However, other decays can lead to complex decay chains such as 

qN^ q{e+£-)N^ , g udC^ ^ udW+N^ . (206) 

With the assumptions that R-parity is conserved and that the iV° is the LSP, all SUSY decay 
chains must end with the N^, which is stable and very weakly interacting. SUSY production pro- 
cesses at hadron colliders then have unbalanced visible momentum, accompanied by multiple jets 
and, possibility, isolated leptons or W and Z bosons. Momentum balance along the beam direc- 
tion cannot be checked at hadron colliders, because fragments of the initial hadrons exit along the 
beam directions, but an imbalance of transverse momentum will be visible and can be a character- 
istic signature of new physics. SUSY events contain this signature and the general large activity 



characteristic of heavy particle production. A simulated event of this type is shown in Fig. 21 



Figure 22 shows a set of estimates given by Tovey and the ATLAS collaboration of the discovery 
potential for SUSY as a function of the LHC luminosity [71] . The most important backgrounds come 
from processes that are themselves relatively rare Standard Model reactions with heavy particle 
production, 

pp {W, Z, tt) + jets . (207) 
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With some effort, we can experimentally normalize and control these backgrounds and reliably 
discovery SUSY production as a new physics process. In the figure, the contours for 5cr excesses of 
events above these backgrounds for various signatures of SUSY events are plotted as a function of 
the so-called 'mSUGRA' parameters. The SUSY models considered are defined as follows: Assume 
gaugino unification with a universal gaugino mass TO1/2 at the grand unification scale. Assume 
also that all scalar masses, including the Higgs boson mass parameters, are unified at the grand 
unification scale at the value ttiq. Assume that the A parameter is universal at the grand unification 
scale; in the figures, the value A = is used. Fix the value of tan (3 at the weak scale. Then it is 
possible to solve for and B, up to a sign, from the condition that electroweak symmetry is broken 
in such a way as to give the observed value of the Z° mass. (I will describe this calculation in 
Section 6.1.) This gives a 4-parameter subspace of the full 24-dimensional parameter space of the 
CP- and flavor-conserving MSSM, with the parameters 

mo , mi/2 , A , tan/3 , sign(^) . (208) 

This subspace is often used to express the results of phenomenological analyses of supersymmetry. In 
interpreting such results, one should remember that this choice of parameters is used for simplicity 
rather than being motivated by physics. 

The figure shows contours below which the various signatures of supersymmetry significantly 
modify the Standard Model expectations. For clarity, the contours of constant squark and gluino 
mass are also plotted. The left-hand plot shows Tovey's results for the missing transverse momentum 
plus multijets signature at various levels of LHC integrated luminosity. It is remarkable that, in the 
models in which the squark or gluino mass is below 1 TeV, SUSY should be discoverable with a data 
sample equivalent to a small fraction of a year of running. The right-hand plot shows the contours for 
the discovery of a variety of SUSY signals, with up to three leptons plus jets plus missing transverse 
momentum, with roughly one year of data at the initial design luminosity. The signals are, as I 
have described, relatively robust with repect to uncertainties in the Standard Model backgrounds. 
This makes it very likely that, if SUSY is really present in Nature as the explanation of electroweak 
symmetry breaking, we will discover it at the LHC. 

The general characteristics of SUSY events also allow us to estimate the SUSY mass scale in 



a relatively straightforward way. In Fig. 23 I show a correlation pointed out by Hinchliffe and 



collaborators |72| between the lighter of the squark and gluino masses and the variable 

4 

M.ff^^T + Y^ETr (209) 

1 

given by the sum of the transverse momenta of the four highest Et jets together with the value of 
the missing transverse momentum. The correlation applies reasonably well to mSUGRA models. 
In other models with smaller mass gaps between the squarks and the lightest neutralino, this relation 
can break down, but M^ff still measures the mass difference between the squark or gluino and the 
Ni ^3]. Some more sophisticated techniques for determining mass scales in SUSY models from 
global kinematic variables are described in [74| . 

5.3 Measurements of the SUSY Spectrum at the LHC 

So far, I have only discussed the observation of the qualitative features of the SUSY model from 
global measures of the properties of events. Now I would like to give some examples of analyses in 
which specific details of the SUSY spectrum are measured with precision at the LHC. The examples 
that I will discuss involve the decay chain 

q^qN^ , N° ^ N^i+£- , (210) 
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Figure 22: Estimates by the ATLAS collaboration of the observability of various signatures of SUSY at 
the LHC. The plots refer to models with grand unification and universal sfermion and gaugino masses Mo 
and Ml/2- The left-hand plot shows the region of this parameter space in which it is possible to detect the 
signature of missing Et plus multiple jets at various levels of integrated luminosity. The right-hand plot 
shows the region of this parameter space in which it is possible to detect an excess of events with one or 
more leptons in addition to jets and missing Et |71| . 
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Figure 23: Correlation between the value of the observable ( |209[ ) and the lighter of the squark and gluino 
masses, from 1721. 



which is typically seen in models in which the gluino is heavier than the squarks and the LSP is 
gaugino-like. 

The decay of the iVj can proceed by any of the mechanisms: 

7v° ^ 7v°z°* , z°* ^ e+r . (211) 

The last line indicates a virtual Z'^, decaying off-shell. In a model with gaugino unification and 
heavy Higgsinos, A'2 is mainly and A^i is mainly 6". Then these modes are preferred in the order 
listed as long as they are kinematically allowed. If the slepton decay is allowed, this is the dominant 
model. Otherwise, the decay to NiZ^ or other open two-body decays dominate. If no t wo-b ody 
decays are open, the A'2 must decay through three-body processes such as the last line of (211 1. 

The decay to an on-shell Z° is hard to work with 75j, but the other two cases can be explored 
in depth. It is useful to begin with the Dalitz plot associated with the 3-body {Ni,£^ ,£~) system. 
Let 

.0 = ^, = = (212) 

m{N2) m{N2) m{N2) 

where the energies are measured in the rest frame of the N2. The three variables are related by 

xi 2:2 = 2 . (213) 
The three-body decay phase space is given by 

^ f dx+ dx^ ; (214) 
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that is, phase space is flat in the variables ( 212[ ). The basic kinematic identities involving the Dalitz 



plot variables are straightforward to work out, especially if we ignore the masses of the leptons. The 
kinematically allowed region is a wedge of the plane bounded by the curves 

x+ + ^ 1 - {m{Ni)/m{N2)Y 
(1 - x+){l - x^) = (TO(iVi)/m(iV2))2 , (215) 

as shown in Fig. [24j^a) . The invariant masses of two-body combinations are given in terms of the 

^^^^Hl-..), (216) 

I am assuming that the Ni is stable and weakly interacting. In this case, the A^i will not be 
observed in the LHC experiments, and also the frame of the N2 cannot be readily determined. 
The only property of this system that is straightforward to measure is the two-body invariant mass 
m{£^£~). So it is interesting to note that the distribution of this quantity distinguishes the first 



and third cases in (211), in the manner shown in Fig. 24 'b). In the case of a two-body decay 
to an intermediate slepton, the decays populate two lines on the Dalitz plot, leading to a sharp 
discontinuity at the kinematic endpoint. In the case of a three-body decay, the events fill the whole 
Dalitz plot, producing a distribution with a slope at the endpoint. With a good understanding of the 
detector resolution in the dilepton invariant mass, these cases can be distinguished experimentally. 

In the three-body case, the endpoint of the dilepton mass distribution is exactly 

m{N2) - m(iVi) , (217) 

so the observable mass distribution gives a precise measurement of this SUSY mass difference. The 
shape of the spectrum has more information. For example, for heavy slepton masses, the shape 
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Figure 25: Distribution of the dilepton invariant mass in two supersymmetry models with 3-body neutralino 
decays: (a.) a model with gaugino-like neutralinos [72], (b.) a model with Higgsino-like neutralinos [73]. In 
the second figure, the dashed curve indicates the m{£^£~) spectrum expected for gaugino-like neutralinos 
with the same mass splitting. 
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Figure 26: Reconstruction of a squark in the model of Fig. 25 a) by combining a dilepton pair at the 
endpoint of the m{£^£~) distribution, the Ni in the same frame with mass determined from kinematics, 
and a &-tagged quark jet. 
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is distinctly different for gaugino-like or Higgsino-like neutralinos. Figure 25 a) shows the dilepton 
mass distribution for an mSUGRA parameter set for which the Ughtest two neutralinos are gaugino- 
like ^72, . Figure [25|^b) shows this distribution for a parameter set in which the two lightest neutralinos 
are Higgsino-like [73]. 

At the endpoint, the dilepton mass is maximal, and this requires that both the dilepton pair and 
the A^i are at rest in the frame of the A^2- By measuring the four- vectors of the leptons, we would 
then know the Ni and N2 four-vectors, up to knowledge of the Ni mass. It is possible to obtain 
this mass approximately from other measurements, for example, from the kinematics of q decays 
directly to A'^i. With this information, we could determine the N2 four- vector. Now the problem of 
missing momentum is solved. By adding observed jets to the A'2 four- vector, it is possible to find 



squarks as resonances [72] . Figure 26 shows the result of such an analysis for the SUSY parameter 



set of Fig. 25 The peak just below 300 GeV is a reconstructed b squark. 



The two-body case of N2 decay is even nicer. In this case, we can see from the right-hand figure 



in Fig. 24 ^b) that the endpoint of the dilepton mass distribution is not located at the mass difference 



(2171 but instead at the smaller value 



m{en-) = m{N2)Jl - Jl - ™^ . (218) 



Figure 27 shows an example of the dilepton spectrum from a SUSY parameter point in this region [63] 
The decay q — > qN2 is also a two-body decay, and there are similar kinematic relations for the upper 
and lower endpoints of the {qtj and {qU) invariant mass distributions. These endpoints are likely 
to be visible in the collider data. Figure [28] shows two jet-lepton mass distributions from a similar 
analysis presented in [76 . In that analysis, it was possible to identify five well-measured kinematic 
endpoints, from which it was possible to solve (in an overdetermined way) for the four masses m{Ni), 
m{t), m{N2), m{q). 

There is one more case of an N2 Ni decay that should be mentioned. If two-body decays of N2 
to sleptons are not kinematically allowed but the decay to Nih° is permitted, this decay to a Higgs 
boson will be the dominant N2 decay. In this case, supersymmetry can provide a copious source of 



Higgs bosons. Figure 29 shows an analysis of a SUSY model in this parameter region [57]. Events 
with multijets and missing transverse energy are selected. In this sample, the mass distribution of 
two 5-quark-tagged jets is shown. The signature of SUSY selects a sample of events in which the 
Higgs boson is visible in its dominant decay to bb. 

There is much more to say about the measurement of SUSY parameters at the LHC. Some more 
sophisticated sets of variables are introduced and applied in [76177] . The question of measuring 
the spins of superparticles is discussed in [7817918 8IJ. And, we have not touched on alternative 
possibilities for the realization of SUSY, with R-parity violation or charged superparticles that are 
observed in the LHC experiments as stable particles. A broader overview of SUSY phenomenology 
at the LHC can be found in the references cited at the beginning of this section. 



6 Electroweak Symmetry Breaking and Dark Matter in the MSSM 

6.1 Electroweak Symmetry Breaking in the MSSM 

In Section 1.2, I motivated the introduction of SUSY with the claim that SUSY could give 
an explanation of electroweak symmetry breaking, and for the presence of weakly interacting dark 
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m„ (GeV) m„ (GeV) 

Figure 27: Dilepton mass distribution in a model with two-body N2 decays, from [63]. The left-hand 
plot shows the dilepton mass distributions for opposite-sign same-flavor dileptons (solid) and for opposite- 
sign opposite-flavor dileptons (dashed). The lower histograms give the estimates of the Standard Model 
background. The right-hand plot shows the difference of the two distributions. 
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Figure 28: Distributions of mass combinations of leptons and high-pT jets showing kinematic endpoints in 
the analysis of [TO|: (a.) the higher m{q£) combination; (b.) the Tn{q£'^£^) distribution. 
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Figure 29: The dijet mass distribution for 2 &-tagged jets at a point in the SUSY parameter space where 
the decay A'^^ — > h'^N^ is dominant, from |67) . 



matter in the universe. Now that we have a detailed understanding of the structure of the MSSM, 
it is time to come back and discuss these issues. 

To present the mechanism of electroweak symmetry breaking in the MSSM, I need to add a term 



to one of the equations that I derived in Section 4.3. In ( 190 ) , I presented the RG equation for the soft 



SUSY breaking scalar mass parameters, including renormalization effects from gauge interactions. 
I remarked that the contributions to this equation from Higgs Yukawa couplings are small for the 
scalars of the first and second generations. However, for the scalars of the third generation, these 
corrections can plan an important role. 



The _F-term interaction 



\ytHu ■ t\ 



(219) 



leads to a contribution to the RG equations for Mt, the mass parameter of t, proportional to Mfj^ 
from the diagram shown in Fig. [30] The value of the diagram is 



(27r)4 B 



(47r)' 



(220) 



A scalar self-energy diagram is interpreted as —i^m}, so this is a negative contribution to . Each 
of the scalar fields {Hu,t,t) gives a similar contribution that renomalizes the soft mass parameter of 
each of the others. For each correction, there is a counting factor from the number of color or SU{2) 
degrees of freedom that run around the loop. There is also a correction to each of the scalar masses 
from the top quark A term. We must also remember that all of these terms add to the positive 
mass correction from the gaugino loops in Fig. |11| of which the gluino loop correction is the most 
important. 



Taking all of these effects into account, we find for the RG equations of the soft mass parameters 
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Figure 30: Diagram contributing a term to the renormalization group equation for the soft mass parameter 
of t proportional to the soft mass parameter for Hu ■ 
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3 • [Af 2 ^ ^^2 ^ ^.^^^ + Al]ml + --- (221) 

The structure is very interesting. The three scalar fields H^, t, and t all receive negative corrections 
to their mass terms as these equations are integrated in the direction of decreasing log Q. If any 
of these mass terms were to become negative, the corresponding field would have an instability 
to develop a vacuum expectation value, and the symmetry of the MSSM would be spontaneously 
broken. The symmetry-breaking we want is that associated with (Hu) =/= 0. However, it seems 

equally possible that we could generate (t) 0, which would break color SU{3), or (<) ^ 0, which 
would break both SU{2) and SU{3). 

If the three mass parameters have similar values at a high mass scale, they race toward negative 



values according to (221 1. But wins the race, and so the theory predicts the symmetry breaking 
pattern that is the one observed. In this way, the MSSM leads naturally to electroweak symmetry 
breaking and realizes the idea that electroweak symmetry breaking is connected to the large value 
of the top quark- Higgs coupling. 



6.2 Higgs Boson Masses in the MSSM 

Once we expect that Af2 < at the weak scale, we can work out the details of the Higgs boson 
spectrum. First, we should write the potential for the Higgs fields Hu, Hd- As in the discussion 
of Sections 4.1 and 4.2, a number of terms need to be collected from the various pieces of the 
Lagrangian. The F terms contriubute 

Vp = ^?{H':*Hl + HTHl) (222) 

The D terms contribute 

Vd = '^^{K*Hl HTH^Y (223) 
The soft SUSY breaking terms contribute 

Vsoft = MluHl*Hl + Mf,dHTH"d - {B^iH'uH'd + h.c.) (224) 
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The sum of these terms gives the complete tree-level Higgs potential. Differentiating this potential 
with respect to H'^ and H^, we obtain the equations that determine the Higgs field vacuum expecta- 
tion values. If we write these equations with the parametrization of the vacuum expectation values 



given in (169), we find 



+ Af|r„ = Bfi cot (3+ -m% cos 2/3 



+ Mffj = Bfi tan /3 -r7i| cos 2^ , 



where m% = {g^ + g''^)v^/4. This system of equations can be solved for id to give 



Af^,-tan^/3A/^, 1 , 



tan^ (3-1 



(225) 



(226) 



This is, for example, the way that we would determine /i in the mSUGRA parameter space described 
in Section 5.2. 

It is interesting to turn this equation around and write it as an equation for mz in terms of the 
SUSY parameters, 

2 _^Mjj^-UT?f3Ml^ 



(227) 



^ tan^ f} - 1 

From this equation, a small value of mz would require a cancellation between the Higgs soft mass 
parameters and ii. The parameter /x sets the mass scale of the Higgsinos, and the Higgs soft mass 
parameters might be related to other masses of the SUSY scalar particles. Thus, if the masses of 
the charginos and neutralinos and, perhaps also, the sleptons are not close to mz, that disparity 
must be associated with an apparently unnatural cancellation between different SUSY parameters. 



If we prohibit a delicate cancellation in (|227|) 

masses 



we put an upper bound on the SUSY partner 
To avoid cancellations in more than two decimal places, must be less than 700 GeV. 
Similarly, we find bounds on the Higgs soft masses, and on the parameters that contribute to these 
masses through the RG equation. This consideration turns out to give a constraint on the gluino 
mass, TO3 < 800 GeV. Assuming gaugino universality, this becomes a condition m2 < 250 GeV 
that restricts the chargino and neutralino masses. A variety of similar naturalness arguments that 
constrain the SUSY scale can be found in |82I83I84] . Though the logic is that of an estimate rather 
than a rigorous bound, this analysis strongly supports the idea that SUSY partners should be hght 
enough to be discovered at the LHC and at the ILC. 

Once we have the Higgs potential and the conditions for the Higgs vacuum expectation values, 
we can work out the masses of the Higgs bosons by expanding the potential around its minimum. 
A first step is to identify the combinations of Higgs fields that correspond to physical Higgs bosons. 
Look first at the charged Higgs bosons. There are two charged Higgs fields in the multiplets Hu, 
Hd- One linear combination of these fields is the Goldstone boson that is eaten by the W boson 
as it obtains mass through the Higgs mechanism. The orthogonal linear combination is a physical 
charged scalar field. If we decompose 



H+ = cos/3i?+ -t-sin/3G+ 
HJ = sin /3H^ + sin /3G^ 



(228) 



where H — (iJ+)*, G — (G"*")*, and (3 is precisely the mixing angle in (1691, it can be seen that 
are the Goldstone bosons and H^ are the physical scalar states. 
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A similar analysis applies to the neutral components of i?^ and H'^ . These are complex- valued 
fields. It is appropriate to decomposed them as 

i7" = ^ (w sin f3 + sin aH° + cos ah° + i cos fiA^ + i sin /JG") 
v2 

iJ^ = ^ (w cos /3 + cos aiJ"- sin a/i° + i sin /3A°-i cos /?G°) (229) 

The components , h° are even under CP; the fields , G" are odd under CP. The componet G" 
is the Goldstone boson eaten by the Z°. The other three fields create physical scalar particles. 

Having identified these fields, we can compute their masses. The formulae for the Higgs masses 
take an especially simple form when they are expressed in terms of the mass of the A^ . For the 
charged Higgs boson 

m|^+ = m\ + . (230) 
For the CP-even scalars, one finds a mass matrix 



m^jsin j3~\-m\cos (3 — (m^ -t- m^) sin /3 cos /3 
— (to^ -I- m|) sin /3 cos cos^ /3 -I- m| sin^ /3 



(231) 



The physical scalar masses m\ and m?^ are the eigenvalues of this matrix, defined in such a way 
that TO^ < ra^jj- The angle a in (229) is the mixing angle that defines these eigenstates. 



Taking the trace of (231 1, we find the relation 

m\ + rn^j = m\ + . (232) 

We can also obtain an upper bound on the lighter Higgs mass m\ by taking the matrix element of 
(2311 in the state (cos sin /?). The bound is a very strong one: 

ml < m| cos^ /? < m| . (233) 

This seems inconsistent with lower bounds on the Higgs boson mass from LEP 2, which exclude 
nih < 114 GeV for the Standard Model Higgs and for most scenarios of SUSY Higgs bosons f8F[^ 
However, the one-loop corrections to the tree-level result (231) give a significant positive correction 

5ml = - ^ sin4 (3 log . (234) 

nmfy mi 

This correction can move the mass of the hP up to about 130 GeV. The detailed summary of the 
radiative corrections to the /i° mass in the MSSM is presented in [88]. A very clear and useful 
accounting of the major corrections can be found in [89 . 

It is possible to raise the mass of the by going outside the MSSM and adding additional SU{2) 
singlet superfields to the model. However, this strategy is limited by a general constraint coming 
from grand unification. The requirement that the Higgs couplings do not become strong up to the 
grand unification scale Hmit the mass of the Higgs to about 200 GeV [90]. It is possible to raise 
the mass of the Higgs further only by enlarging the Standard Model gauge group or adding new 
thresholds that affect unification |91l92j . 

In the MSSM, we can easily have the situation in which mA S> mh- In this limit, the couplings 
of the are very close to those of the Standard Model Higgs boson, and the A^, and are 
almost degenerate. If tan /3 ^ 1, the heavy neutral Higgs bosons decay dominantly to bb and t~^t~. 

Much more about the phenomenology of Higgs bosons in supersymmetry can be found in [93194] . 



5 Some exceptional Higgs decay schemes that escape these bounds are considered in | 86I87| . 
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6.3 WIMP Model of Dark Matter 



Now we turn to the second problem highlighted in the Introduction, the problem of dark matter in 
the universe. It has been known from many astrophysical measurements that the universe contains 
enormous amounts of invisible, weakly interacting matter. For an excellent review of the classic 
astrophysical evidence for this dark matter, see |95j . 

In the past few years, measurements of the cosmic microwave background have given a new 
source of evidence for dark matter. Since this data comes from an era in the early universe before 
the formation of any structure, it argues strongly that the invisible matter is not made of rocks or 
brown dwarfs but is actually a new, very weakly interacting form of matter. These measurements 
also determine quite accurately the overall amount of conventional and dark matter in the universe. 
Let pb, pn, and pA be the large-scale energy densities of the universe from baryons, dark matter, and 
the energy of the vacuum. The data from the microwave background tells us that pb + pn + Pa = Pc, 
the 'closure density' corresponding in general relativity to a flat universe, to about 1% accuracy. If 

— Pi/Pc, the most recent data from the WMAP experiment and other sources gives |96l97j 

r^b = 0.042 ± 0.003 r^AT = 0.20 ± 0.02 r^A = 0.74 ± 0.02 . (235) 

These results present a double mystery. We do not know what particle the dark matter is made of, 
and we do not have any theory that explains the observed magnitude of the vacuum energy or 'dark 
energy'. 

I believe that supersymmetry will eventually play an essential role in solving the problem of dark 
energy. In ordinary quantum field theory, the value of the vacuum energy is quartically divergent, 
so the problem of computing the vacuum energy is not even well-posed. In supersymmetry, there 
is at least a well-defined zero of the energy associated with exact supersymmetry, which implies 
(0| H |0) — 0. Unfortunately, in most of today's models of supersymmetry, the vacuum energy is 
set by the SUSY breaking scale. This gives A ~ (10^^ GeV)^, about 80 orders of magnitude larger 
than the observed value of the vacuum energy. From this starting point, A must be fine-tuned to 
the scale of eV^. This is an important problem that needs new insights which, however, I will not 
provide here. 

On the other hand, supersymmetry offers a very definite solution to the problem of the origin 
of dark matter. We have already noted in Section 3.4 that it is straightforward to arrange that 
the lightest supersymmetric particle can be absolutely stable. If this particle were produced in the 
early universe, some density of this type of matter should still be present. In most, but not all, 
regions of parameter space, the lightest supersymmetric particle is neutral. Candidates include the 
lightest neutralino, the lightest sneutrino, and the gravitino. In the remainder of these lectures, 
I will concentrate on the case in which the lightest neutralino is the dark matter particle. For a 
discussion of the other candidates, see [98] . 

To begin our discussion, I would like to estimate the cosmic density of dark matter in a more 
general context. Let me make the following minimal assumptions about the nature of dark matter, 
that the dark matter particle is stable, neutral, and weakly interacting. To these properties, I 
would like to add one more, that dark matter particles can be created in pairs at sufficiently high 
temperature, and that, at some time in the early universe, dark matter particles were in thermal 
equilibrium. I will refer to a particle satisfying these assumptions as a 'weakly interacting massive 
particle' or WIMP. The assumption of thermal equilibrium is a strong one that is not satisfied even 
in many models of supersymmetric dark matter. For some exceptions, see |99ll00j . However, let us 
see what implications follow from these assumptions. 

The assumption that WIMPs were once in thermal equilibrium provides a definite initial condition 
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from which to compute the current density of dark matter. In thermal equihbrium at temperture 
T, we have for the number density of dark matter particles 



-^. = J^.^-^Tr'<^-'- . (236) 

where g is the number of spin degrees of freedom of the massive particle. As the universe expands, 
the temperature of the universe deccreases and the rate of WIMP pair production becomes very 
small. But the rate of dark matter pair annihilation also becomes small as the WIMPs separate 
from one another. 

The expansion of the universe is governed by the Hubble constant H = a/a, where a is the scale 
factor. Einstein's equations imply that 

H'='i^- (237) 

In a radiation-dominated universe where is the number of relativistic degrees of freedom, p = 
7T^g^,T'^/30. Then H is proportional to T^. In a radiation-dominated universe, the temperature 
red-shifts as the universe expands, so that T ^ a^^. Combining this relation with the equation 
H = a/a ^ T^, we find t ^ ~ a^, that is, a ^ t^^^ or a/a — l/2t. Setting this expression 



equal to the explict form of H in (237), we find a detailed formula for the time since the start of the 



radiation-dominated era for cooling to a temperature T, 

-(i^)"">^ 

The evolution of the WIMP density is described by the Boltzmann equation 

^ = -3Hn - {av) [n" - n\) , (239) 

where H is the Hubble constant, a is the NN annihilation cross section — which appears thermally 
averaged with the relative velocity of colliding WIMPs — and n^n is the equilibrium WIMP density 



(2361. Assume, just for the sake of argument, that the temperature T is of the order of 100 GeV. 
At this temperature, the Hubble constant has the magnitude H ~ 10~^^r, so the expansion of the 
universe is very slow on the scale of typical elementary particle reactions. However, when T becomes 
less than the WIMP mass m, the WIMP density is exponentially suppressed and so the collision 
term in the Boltzmann equation is also very small. These two terms are of the same size at the 
freezeout temperature Tp satisfying 

e-™/^^ . (240) 

nipim [av) 

At temperatures below Tp, we may neglect the production of WIMPs in particle collisions. The 
WIMP density is then determined by the expansion of the universe and the residual rate of WIMP 
pair annihilation. Maybe it is more appropriate to think of Tp as the temperature at which a WIMP 
density is frozen in. To determine the freezeout temperature, we take the logarithm of the right-hand 



side of ( 240 1 . The result depends only on the order of magnitude of the annihilation cross section. 



For any interaction of electroweak strength, 

^P ^Tp/m^ 1/25 . (241) 

This physical picture suggests a way to estimate the cosmic density of WIMP dark matter. We 
can take as our initial condition the thermal density of dark matter at freezeout. We then integrate 
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the Boltzmann equation, ignoring the term proportional to associated with the production of 
WIMP pairs [TO^ . 

In analyzing the Boltzmann equation, it is useful normalize the particle density n of dark matter 
to the density of entropy s. Since the universe expands very slowly, this expansion is very close to 
adiabatic. Then entropy is conserved, 

^ = -3Hs . (242) 
In a radiation-dominated universe, s = 2Tr^g^T^ /45. Now define 



Y = 



T 



(243) 



the latter as in (2411. Using the expression (2381, we can convert the evolution in time to an 



evolution in temperature or in ^. Applying these changes of variables and dropping the nj^ term 



the Boltmann equation ( 239 1 rearranges to the form 



dY 



C (av) Y^ 



where 



C 



V 45 / 



1/2 



mmpi . 



(244) 



(245) 



Let Yp be the value of y at ^ = S^p. If we assume that (av) is approximately constant, since we 
are at temperatures close to threshold, it is straightforward to integrate this equation to ^ = 0, 
corresponding to late times. 



Y- 



Y^ 



(246) 



The second term typically dominates the first. Then we can put back the value of C in (2451 and 



write the final answer in terms of the ratio of the mass density of dark matter to the closure density 
^Im ~ nmi^ / pc- In this way, we find 



Q.N — 



So ( 45 



Pc \T^g 



1/2 



1 



1 



(247) 



where so is the current entropy density of the universe. Turner and Scherrer observed that this 
formula gives a value of fi^v that is usually within 10% of the result from exact integration of the 
Boltzmann equation |102j . If ( sigmav) has a significant dependence on temperature, the derivation 
is still correct with the replacement 



CM 



d^{av) (0 



(248) 



in the denominator of the last term in ( 247 1 



This is a remarkable relation. Almost every factor in this relation is known from astrophysical 



measurements. The left-hand side is given by (2351. On the right-hand side, the entropy density of 



the universe is dominated by the entropy of the microwave background photons and can be computed 
from the microwave background temperature. The closure density is known from the measurement 
of the Hubble constant and the observation that the universe is flat. The parameters f/* and (,f 
are relatively insensitive to the strength of the annihilation cross section, with values ^ 100, 
£,p ^ 1/25. The mass of the WIMP does not appear explicitly in (2471. We can then solve for (av). 
The result is 

(av) = 1 pb . (249) 
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This is the value of a typical electroweak cross section at energies of a few hundred GeV. If we 
convert this value to a mass M of an exchanged particle using the formula 

(250) 



MP 



the value (249) corresponds to M = 100 GeV. 



I consider this a truly remarkable result. From a purely astrophysical argument, relying on quite 
weak and general assumptions, we arrive at the conclusion that there must be new particles at the 
hundred GeV energy scale. It is probably not a concidence that this argument leads us back to the 
mass scale of electroweak symmetry breaking. 

In our study of supersymmetry, we have found an argument from the physics of electroweak 
symmetry breaking that predicts the existence of dark matter. As I discussed at the beginning of 
these lectures, models that explain electroweak symmetry breaking are complex. They typically 
involve many new particles. It is easily arranged that the lightest of the new particles is neutral. 
In supersymmetry, there is a reason why the new particles are likely to carry a conserved quantum 



number (149 1. Other models of electroweak symmetry breaking, such as the extra dimensional 
and little Higgs models discussed in Section 1.2, have their own reasons to have a complex particle 
spectrum and discrete symmetries. Then these models lead in their own ways to WIMPs at the 
hundred GeV mass scale. 

A slight extension of this argument adds more interest. In supersymmetry, the sector of new 
particles includes particles with QCD color. Since the top quark probably plays an essential role in 
the mechanism of electroweak symmetry breaking, it is very likely that, in any model, some of the 
new particles will carry color. If these particles have masses below 1 TeV, they have large (10 pb) 
pair-production cross sections at the LHC. These particles will then decay to the dark matter particle, 
producting complex events with several hard jets and missing transverse momentum. These mild 
assumptions thus lead to the conclusion, from any model that follows this general line of argument, 
that we should expect exotic events with multiple jets and missing transverse momentum to appear 
with pb cross sections at the LHC. 



6.4 Dark Matter Annihilation in the MSSM 

This argument of the previous section gives a very optimistic conclusion for the discovery of new 
physics at the LHC. However, we have already discussed that the first observation of supersymmetry 
or another model of new physics will only be the first step in a lengthy experimental program. Once 
we know that superparticles or other new particles exist, we will need to study them in detail to 
learn their detailed interactions and, eventually, to work out the underlying Lagrangian that governs 
their behavior. As we have already discussed in Section 3.5 and 4.3, this Lagrangian can give us a 
clue to the nature of the ultimate theory at very short distances. 

The study of dark matter intersects this program in an interesting way. In principle, once we 
have discovered supersymmetric particles, we can try to measure their properties and see if these 
coincide with the properties required from astrophysical detections of dark matter. As we have seen 
in Section 5.3, the LHC experiments expect to measure the mass of the LSP to about 10% accuracy. 
These measurements can hopefully be compared to mass measurements at the 20% level that can 
be expected from astrophysical dark matter detection experiments [1031104] . We would also wish 
to find out whether the annihilation cross section (av) that is predicted from the supersymmetry 



parameters measured at colliders agrees with the value ( 249 ) required to predict the observed WIMP 
relic density. This comparison turns out to depend in a complex way on the parameters of the 
underlying supersymmetry theory. 
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Figure 32: Two possible spin configurations for neutralino annihilation: (a.) spin 0; (b.) spin 1. Because of 
Fermi statistics, the latter state does not exist in the S-wave. 



To begin our discussion of the annihilation cross section, we can make a simple model of neutralino 
annihilation and see how well it works. We have seen in Section 4.3 that the right-handed sleptons 
are often the lightest charged particles in the supersymmetry spectrum. Consider, then, an idealized 
parameter set in which the neutralino is a pure bino and pair annihilation is dominated by the 
slepton exchange diagrams shown in Fig. |31[ (Away from the pure bino case, there are also s- 
channel diagrams with , hP , H'^, A*^.) In this special limit, the annihilation cross section is given 

by 

da 9 



dcos0 



1 


2 


1 1 






m~ — t mi. — u 



(251) 



where tojv is the iVi mass. The relative velocity v appears due to the flux factor in the cross section; 
this factor cancels in av. I have ignored the lepton masses. This expression is of the order of 
(2501 with M ~ m^v, except for one unfortunate feature: At threshold, t = u and the cross section 



vanishes. This leads to a severe suppression, by a factor of 



'TV 



?7iH,+ m^ 



N 



(252) 



which is at least of order So the relic density estimated in this simple way is too large by 

about a factor of 10. 

There is an interesting physics explanation for the vanishing of this cross section at thresh- 
old |105j . Neutralinos are spin-| fermions, and we might guess from this that, near threshold, they 
would annihilate in the S-wave either in a spin or in a spin 1 state. The two spin configurations 



are shown in Fig. 32 However, because the neutralino is a Majorana fermion and therefore its own 
antiparticle, an S-wave state of two neutralinos must be antisymmetric in spin. Hence, the spin 1 
S-wave state does not exist However, as we know from pion decay, a spin state can convert to 
a pair of light leptons only with a helicity flip. Thus, there is an annihilation cross section from 
the spin S-wave only when lepton masses are included, and even then with the suppression factor 
m'j/mff, which is 10^'* even for r+r^ flnal states. 
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Figure 33: Three mechanisms for obtaining a sufficiently large annihilation cross section to give the observed 
density of neutralino dark matter: (a.) gaugino-Ifiggsino mixing, opening the annihilation channels to 
W^W~ and Z'^ Z° , (b.) resonance annihilation through the Higgs boson A'^ , (c.) co-annihilation with 
another supersymmetric particle, here taken to be a I. 



To obtain a realistic value for the neutralino relic density, we have to bring in more complicated 
mechanisms of neutralino annihilation. These mechanisms are not difficult to find in various regions 
of the large supersymmetry paramet er space PD6 107 108^ . We need to look for annihilation pro- 
cesses that can proceed in the S-wave with full strength. Three possible mechanisms are shown in 

Fig.m 

Pairs of neutralinos can annihilate in the S-wave into vector bosons. The bino does not couple 
to W or Z pairs, but if the lightest neutralino has Higgsino or wino content, this reaction can be 
important. For charginos of mass about 200 GeV, this annihilation cross section can be 50 pb for a 
pure wino or Higgsino, so only a modest content of these states is needed to give a cross section of 
1 pb. 

The s-channel exchange of a Higgs boson can provide a mechanism for neutralino annihilation in 
the spin S-wave. Because this state is CP-odd, it is the boson that is relevant here. If is close 
to the neutralino threshold 2m n, the cross section has a resonant enhancement. Note that the A^i 
annihilation vertex to A arises as a Higgs-Higgsino-gaugino Yukawa term, so this vertex is nonzero 
only if A^i has both gaugino and Higgsino content. If = 2tojv, the resonance enhancement is at 
full strength and the cross section can be as large as 50 pb. Thus, it is A boson masses about 20 



GeV above or below the threshold that give the desired cross section ( 249 1 



The final mechanism shown in the figure is coannihilation. As we have discussed, the freezeout 
of the A^i occurs at a temperature given by T/m^ ^ 1/25. So if there is another particle in the 
supersymmetry spectrum that is within 4% of the A^i mas^, this state will have a number density that 
remains in equilibrium with the number density of the A^i . If this particle has S-wave annihilation 
reactions, those reactions can be the dominant mechanisms for the annihilation of supersymmetric 
particles. For a light slepton, the reactions 

J- + N^^r + , r +1- ^ r + r (253) 

can give significant S-wave annihilation. In |106ll09j . the lighter stau is invoked as the coannihilating 
particle. In |110j . the lighter top squark is invoked as the coannihilating state. If the lightest 
neutralinos and charginos are Higgsino- like, chargino coannihilation can also be important. 

It is, then, a complex matter to predict the neutralino relic density from microscopic physics. 
We will first need to learn what particles in the supersymmetry spectrum play the dominant role 
as particle exchanged in annihilation reactions or as coannihilating species. We will then need to 
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Figure 34: Scatter plot of SUSY parameter points consistent with data from the LHC in the analysis of 
the parameter set LCC2 from The horizontal axis show the value of 57]v at each parameter point. 

The vertical axes show polarized-beam cross sections measurable at the ILC, in fb: (a.) CT(e^eJ — ^ C^Ci), 
(b.) cr(e^ej N2Nli). The colored bands show the ±la region allowed after the ILC cross section 
measurements. 



measure the couplings and mixing angles of the important particles, since the dominant annihilation 
diagrams depend sensitively on these. 

Some examples of how measurements at the LHC and ILC can accumulate the relevant infor- 



mation are described in Figure 34 shows a part of the analysis of this paper for a particular 

SUSY model in which the dominant annihilation reactions are NiNi W~^W~ , Z^Z^. As a first 
step, the authors constructed numerous supersymmetry parameter sets that were consistent with the 
mass spectrum of this model as it would be measured at the LHC. These parameter sets included a 
variety of models in which the LSP was dominantly bino and wino. The figure shows scatter plots of 
the predictions of these models with ILC cross sections for neutralino and chargino pair production 
on the vertical axis and fi^r on the horizontal axis. The two cross sections clearly separate the bino- 
and wino-like solutions. The second of these cross sections is the polarized reaction of chargino 
pair production for which the cross section is displayed in Fig. |19[ The horizontal lines represent 
the accuracy of the measurements of these cross sections expected at the ILC. These measurements 
select the bino solution and also play an important role in fixing the bino-Higgsino mixing angle 



which is a crucial input to the annihilation cross sections. In Fig. 35 I show the distribution of 
predictions for ilpf expected for this model, in the analysis of from the data on SUSY particles 
that would be obtained from the LHC, from the ILC at a center-of-mass energy of 500 GeV, and 
from the ILC at a center-of-mass energy of 1000 GeV. 



The similar summary plot for another of the models considered in [111] is shown in Fig. 36 
The model considered in this analysis is one in which the neutralino relic density is set by stau 
coannihilation. In this model, the stau would be discovered at the LHC, and the stau-neutralino 
mass difference would be measured to about 10% accuracy at the 500 GeV ILC. However, the 
annihilation reactions also depend on mixing angles and on the value of tan /?. In this scenario, these 
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Figure 35: Summary plot for the prediction of f2]v from collider data for the SUSY parameter set LCC2 
considered in The three curves show the likelihood distributions for the prediction of f2jv using data 

from the LHC, the ILC at 500 GeV, and the ILC at 1000 GeV. 



are determined only by ILC measurements of some of the heavier states of the SUSY spectrum. 

Collider measurements of the SUSY spectrum can also be used to constrain cross sections of 
the WIMP that are important for experiments that seek to detect dark matter, for example, the 
neutralino-proton cross section and the cross section for neutralino pair annihilation to gamma rays. 
If we can accurately predict these cross sections from collider data, the information about the SUSY 
spectrum that we learn from colliders will feed back into the astrophysics of dark matter. Some 
numerical examples that illustrate this are presented in |lllj . 



7 Conclusions 

In these lectures, I have given an overview of supersymmetry and its application to elementary 
particle physics. In the early sections of this review, I presented the formalism of SUSY and explained 
the rules for constructing supersymmetric Lagrangians. Our discussion then became more concrete, 
focusing on the mass spectrum of the MSSM and the properties of the particle states of the MSSM 
spectrum. This led us to a discussion of the experimental probes of this spectrum and the possibility 
of measurement of the parameters of the supersymmetric Lagrangian. 

This possibility is now coming very near. As I have discussed in the last sections of this review, 
supersymmetry gives concrete answers to the major questions about elementary particle physics 
that we expect to be addressed at the hundred GeV scale — the questions of the origin of electroweak 
symmetry breaking and the identity of cosmic dark matter. In the next year, the LHC will begin 
to explore the physics of this mass scale. Supersymmetry is one candidate for what will be found. I 
hope that, after studying these lectures, you will agree that the picture provided by supersymmetry 
is highly plausible and even compelling. 
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Whatever explanations we will learn from the LHC data, our investigation of it will follow the 
general paradigm that I have described here. In successive stages, we will use data from the LHC 
and the ILC to learn the mass spectrum of new particles that are revealed at the LHC, to determine 
their quantum numbers and couplings, and to reconstruct their underlying Lagrangian. On the 
basis of the detailed studies of this program that have been carried out for the MSSM, we have the 
expectation that we will be able to learn the underlying theory of the new particles and to test the 
specific explanations that this theory gives for the mysteries of the fundamental interactions. 

Is supersymmetry just an attractive theory, or is it a part of the true description of elementary 
particles? We are about to find out. 
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